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Abstract 

We investigate the steepening of the magnetic fluctuation power 
law spectra observed in the inner solar wind for frequencies higher 
than 0.5 Hz. This high frequency part of the spectrum may be at- 
tributed to dispersive nonlinear processes. In that context, the long- 
time behavior of weakly interacting waves is examined in the frame- 
work of three-dimensional incompressible Hall magnetohydrodynamic 
(MHD) turbulence. The Hall term added to the standard MHD equa- 
tions makes the Alfven waves dispersive and circularly polarized. We 
introduce the generalized Elsasser variables and, using a complex he- 
licity decomposition, we derive for three-wave interaction processes 
the general wave kinetic equations; they describe the nonlinear dy- 
namics of Alfven, whistler and ion cyclotron wave turbulence in the 
presence of a strong uniform magnetic field Bq&\\. Hall MHD tur- 
bulence is characterized by anisotropies of different strength: (i) for 
wavenumbers kdi ^> 1 (di is the ion inertial length) nonlinear trans- 
fers are essentially in the direction perpendicular (_L) to Bo; (ii) for 
kdi "C 1 nonlinear transfers are exclusively in the perpendicular direc- 
tion; (iii) for kdi ~ 1, a moderate anisotropy is predicted. We show 
that electron and standard MHD turbulence can be seen as two fre- 
quency limits of the present theory but the standard MHD limit is 
singular; additionally, we analyze in detail the ion MHD turbulence 
limit. Exact power law solutions of the master wave kinetic equa- 
tions are given in the small and large scale limits for which we have, 
respectively, the total energy spectra E(k±, k\\) ~ A: ± 5//2 |/c|| |~ 1//2 and 
E(k±,hi) ~ kj 2 . An anisotropic phenomenology is developed to de- 
scribe continuously the different scaling laws of the energy spectrum; 



1 



one predicts E(k±,ku) ~ /cJ 2 |A;|| | _1 / 2 (1 + k\d1)~ l 1 4 . Nonlocal inter- 
actions between Alfven, whistler and ion cyclotron waves are investi- 
gated; a non trivial dynamics exists only when a discrepancy from the 
equipartition between the large scale kinetic and magnetic energies 
happens. 

1 Introduction 

Spacecraft observations of the inner solar wind, i.e. for heliocentric distances 
less than 1AU, show magnetic and velocity fluctuations over a broad range 
of frequencies, from 10~ 5 Hz up to several hundred Hz (see e.g. Coleman, 
1968; Belcher and Davis, 1971; Matthaeus and Goldstein, 1982; Roberts et 
al., 1987; Grappin et al., 1990; Burlaga, 1991; Leamon et al., 1998). These 
fluctuations possess many properties expected of fully developed weakly com- 
pressible magnetohydrodynamic (MHD) turbulence (Goldstein and Roberts, 
1999). For that reason the interplanetary medium is often seen as a vast lab- 
oratory for studying many fundamental questions about turbulent plasmas. 

Compressed interaction regions produced by velocity differences are clearly 
observed in the outer solar wind. In regions where fast streams overtake slow 
streams, the spectrum of the large scale magnetic field fluctuations follows 
a f~ 2 frequency power law (Burlaga et al., 1987) which was shown to be a 
spectral signature of jumps (Roberts et al., 1987). Compressive effects are 
however weaker in the inner solar wind. For example, the normalized den- 
sity fluctuations near the current sheet at 0.3AU are often smaller than 5% 
(Bavassano et al., 1997). This tendency is confirmed by indirect measure- 
ments through radio wave interplanetary scintillation observations at helio- 
centric distances of 16 — 26 solar radius (Spangler, 2002). Waves and turbu- 
lence - the subject of this paper - are better observed in the pure/polar wind 
where generally the density fluctuations are weaker than in the current sheet. 
Therefore, the inner solar wind may be seen mainly as a weakly compressible 
medium. 

The turbulent state of the solar wind was suggested by [21] who reported 
a power law behavior for energy spectra with spectral indices lying between 
— 1 and —2. The original signals being measured in time, these observa- 
tional spectra are measured in frequency. Since the solar wind is super- 
sonic and super-Alfvenic, the Taylor "frozen-in flow" hypothesis is usually 
used to connect directly a frequency to a wavenumber which allows even- 
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tually a comparison with theoretical predictions. Note, however, that for 
anisotropic turbulence a relevant comparison with theoretical predictions, 
like those made in this paper, are only possible if a three dimensional en- 
ergy spectrum is accessible by in situ measurements, a situation that is not 
currently achieved (see, however, Matthaeus et al., 2005; see also the last 
Section). More precise measurements (see e.g. Matthaeus and Goldstein, 
1982) revealed that the spectral index at low frequency is often about — 1.7 
which is closer to the Kolmogorov prediction (Kolmogorov, 1941) for neutral 
fluids (—5/3) rather than the Iroshnikov-Kraichnan prediction (Iroshnikov, 
1963; Kraichnan, 1965) for magnetized fluids (—3/2). Both predictions are 
built, in particular, on the hypothesis of isotropic turbulence. However, the 
presence of Alfven waves in the fast solar wind attests that the magnetic 
field has a preferential direction which is likely at the origin of anisotropic 
turbulence provided that the amount of counter-propagating Alfven waves 
is enough. Indeed, in situ measurements of cross helicity show clearly that 
outward propagative Alfven waves are the main component of the fast solar 
wind at short radial distances (Belcher and Davis, 1971) but they become 
less dominant beyond 1AU. Since pure Alfven waves are exact solutions of 
the ideal incompressible MHD equations (see e.g. Pouquet, 1993), nonlinear 
interactions should be suppressed if only one type of waves is present. The 
variance analysis of the magnetic field components and magnitude shows 
clearly that the magnetic field vector of the fast solar wind has a varying di- 
rection but with only a weak variation in magnitude (see e.g. Forsyth et al., 
1996ab). Typical values give a ratio of the normalized variance of the field 
magnitude smaller than 10% whereas for the components it can be as large 
as 50%. In these respects, the interplanetary magnetic field may be seen as 
a vector liying approximately around the Parker spiral direction with only 
weak magnitude variations (Barnes, 1981). Solar wind anisotropy with more 
power perpendicular to the mean magnetic field than that parallel is pointed 
out by data analyses (Belcher and Davis, 1971; Klein et al., 1993) with a 
ratio of power up to 30. From single-point spacecraft measurements it is 
however not possible to specify the exact three-dimensional form of the spec- 
tral tensor of the magnetic or velocity fluctuations. In absence of such data, 
[15J with a quasi two-dimensional model, in which wave vectors are nearly 
perpendicular to the large-scale magnetic field, argued that about 85% of so- 
lar wind turbulence possesses a dominant 2D component. Additionally, solar 
wind anisotropies is detected through radio wave scintillations which reveal 
that density spectra close to the Sun are highly anisotropic with irregularities 
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stretched out mainly along the radial direction (Armstrong et al., 1990). 

Most of the papers dealing with interplanetary turbulence tend to focus 
on a frequency inertial range where the MHD approximation is well satis- 
fied. It is the domain where the power spectral index is often found to be 
around the Kolmogorov index. During the last decades several properties 
of the solar wind turbulence have been understood in this framework. Less 
understood is what happens outside of this range of frequencies. At lower 
frequencies (< 10~ 5 Hz) flatter power laws are found in particular for the fast 
solar wind with indices close to —1 and even less (in absolute value) for the 
smallest frequencies. These large scales are often interpreted as the energy 
containing scales. But the precise role of the low solar corona in the gen- 
eration of such scales and the origin of the evolution of the spectral index 
when the distance from the Sun increases are still not well understood (see 
e.g. Velli et al., 1989; Horbury, 1999). For frequencies higher than 0.5Hz a 
steepening of the magnetic fluctuation power law spectra is observed over 
more than two decades (Coroniti et al., 1982; Denskat et al., 1983, Leamon 
et al., 1998) with a spectral index on average around —3. Note that the 
latest analysis made with the Cluster spacecraft data reveals a less steep 
spectral index about —2.12 (Bale et al., 2005). This new range, exhibiting 
a power law, is characterized by a bias of the polarization suggesting that 
these fluctuations are likely to be right-hand polarized, outward propagating 
waves (Goldstein et al., 1994). Various indirect lines of evidence indicate 
that these waves propagate at large angles to the background magnetic field 
and that the power in fluctuations parallel to the background magnetic field 
is much less than the perpendicular one (Coroniti et al., 1982; Leamon et 
al., 1998). For these reasons, it is thought (e.g. Stawicki et al., 2001) that 
Alfven - left circularly polarized - fluctuations are suppressed by proton cy- 
clotron damping and that the high frequency power law spectra are likely 
to consist of whistler waves. This scenario proposed is supported by direct 
numerical simulations of compressible 2|D Hall MHD turbulence (Ghosh et 
al., 1996) where a steepening of the spectra is found and associated with the 
appearance of right circularly polarized fluctuations. It is plausible that what 
has been conventionally thought of as a dissipation range is actually a disper- 
sive or inertial range and that the steeper power law may be due to nonlinear 
wave processes rather than dissipation (see e.g. Krishan and Mahajan, 2004). 
Under this new interpretation, the resistive dissipation range of frequencies 
may be moved to frequencies higher than the electron cyclotron frequency. A 
recent study (Stawicki et al., 2001) suggests that the treatment of the solar 
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wind dispersive range should include magnetosonic/whistler waves since it 
is often observed that the high frequency fluctuations of the magnetic field 
are much smaller than the background magnetic field (see also Forsyth et 
al., 1996ab). It seems therefore that a nonlinear theory built on weak wave 
turbulence may be a useful point of departure for understanding the detailed 
physics of solar wind turbulence. 

It is well known that the presence of a mean magnetic field plays a fun- 
damental role in the behavior of compressible or incompressible MHD tur- 
bulence (Montgomery and Turner, 1981; Shebalin et al. 1983; Oughton et 
al, 1994; Matthaeus et al, 1996; Kinney and McWilliams, 1998; Cho and 
Vishniac, 2000; Oughton and Matthaeus, 2005). The main effect is that the 
mean magnetic field renders the turbulence quasi-bidimensional with a non- 
linear transfer essentially perpendicular to its direction. Note that in such 
a situation, it is known that compressible MHD behaves very similarly to 
reduced MHD (see e.g. Dmitruk et al, 2005). This property is generalized 
to 2|D compressible Hall MHD for high and low beta (3 plasma simulations 
(Ghosh and Goldstein, 1997) for which strong anisotropies are also found 
when a strong mean magnetic field is present. It was suggested that the 
action of the Hall term is to provide additional suppression of energy cas- 
cades along the mean field direction and incompressibility may not be able to 
reproduce the dynamics seen in the simulations. This conclusion is slightly 
different to what we find in the present study where we show that incom- 
pressibility in Hall MHD allows to reproduce several properties observed like 
anisotropic turbulence. Incompressible Hall MHD is often used to under- 
stand, for example, the main impact of the Hall term on flowing plasmas 
(Ohsaki, 2005), or in turbulent dynamo (Mininni et al, 2003). In the present 
paper, we derive a weak wave turbulence formalism for incompressible Hall 
MHD in the presence of a strong external mean magnetic field, where Alfven, 
ion cyclotron and whistler/electron waves are taken into account. One of the 
main results is the derivation of the wave kinetic equations at the lowest 
order, i.e. for three- wave interaction processes. For such a turbulence, it is 
possible to show, in particular, a global tendency towards anisotropy with 
nonlinear transfers preferentially in the direction perpendicular to the exter- 
nal magnetic field. Hall MHD wave turbulence theory describes a wide range 
of frequencies, from the low frequency limit of pure Alfvenic turbulence to 
the high frequency limit of whistler wave turbulence for which the asymp- 
totic theories have been derived recently (Galtier et al., 2000-2002; Galtier 
and Bhattacharjee, 2003). By recovering both theories as two particular lim- 
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its, we recover all the well-known properties associated. Ion cyclotron wave 
turbulence appears as a third particular limit for which we report a detailed 
analysis. The energy spectrum of Hall MHD is characterized by two inertial 
ranges, which are exact solutions of the wave kinetic equations, separated by 
a knee. The position of the knee corresponds to the scale where the Hall term 
becomes sub/dominant. We develop a single anisotropic phenomenology that 
recovers the power law solutions found and makes the link continuously in 
wavenumbers between the two scaling laws. A nonlocal analysis performed 
on the wave kinetic equations reveals that a non trivial dynamics between 
Alfven, whistler and ion cyclotron waves happens only when a discrepancy 
from the equipartition between the large scale kinetic and magnetic energies 
exists. We believe that the description given here may help to better under- 
stand the inner solar wind observations and, in particular, the existence of 
dispersive/non-dispersive inertial ranges. 

The organization of the paper is as follows: in Section 2, we discuss 
about the approximation of incompressible Hall MHD and the existence of 
transverse, circularly polarized waves; we introduce the generalized Elsasser 
variables and the complex helical decomposition. In Section 3, we develop the 
wave turbulence formalism and we derive the wave kinetic equations for three- 
wave interaction processes. Section 4 is devoted to the general properties of 
Hall MHD wave turbulence. Section 5 deals with the small and large scale 
limits respectively of the master equations of wave turbulence. In particular, 
we give a detail study of the ion cyclotron wave turbulence. In Section 6, 
we derive an anisotropic heuristic description which, in particular, makes 
the link between the previous predictions. In Section 7, we analyze nonlocal 
interactions between Alfven, whistler and ion cyclotron waves. Section 8 
is devoted to the possible sources of anisotropy. In Section 9, we discuss 
about the domain of validity of the theory and the difference between wave 
and strong turbulence. Finally, we conclude with a summary and a general 
discussion in the last Section. 

2 Hall magnetohydrodynamics approximation 

2.1 Generalized Ohm's law 

Hall MHD is an extension of the standard MHD where the ion inertia is 
retained in Ohm's law. The generalized Ohm's law, in SI unit, is then given 
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by 

E + VxB ^ =/ior? J, (1) 

ne 

where E is the electric field, V is the plasma flow velocity, B is the magnetic 
field, J is the current density, n is the electron density, e is the magnitude 
of the electron charge, uo is the permeability of free space and rj is the 
magnetic diffusivity. The Hall effect, represented by the last term in the left 
hand side of the generalized Ohm's law, becomes relevant when we intend 
to describe the plasma dynamics up to length scales shorter than the ion 
inertial length di (di = c/u P i, where c is the speed of light and u pi is the 
ion plasma frequency) and time scales of the order, or shorter, than the ion 
cyclotron period oo^ 1 . It is one of the most important manifestations of the 
velocity difference between electrons and ions when kinetic effects are not 
taken into account. The importance of the Hall effect in astrophysics has 
been pointed out to understand, for example, the presence of instabilities in 
protostellar disks (Balbus and Terquem, 2001), the magnetic field evolution 
in neutron star crusts (Goldreich and Reisenegger, 1992; Cumming et al, 
2004), impulsive magnetic reconnection (see e.g. Bhattacharjee, 2004) or the 
formation of filaments (see e.g. Passot and Sulem, 2003; Dreher et al., 2005). 



2.2 Incompressible Hall MHD equations 

The inclusion of the Hall effect in the Ohm's law leads, in the incompressible 
case, to the following Hall MHD equations: 

V-V = 0, (2) 

^+V-VV = -VP»+B-VB + ^V 2 V, (3) 
ot 

— +V-VB = B-VV-diVx [(Vx B) x B] + V V 2 B , (4) 

V-B = 0, (5) 

where B has been normalized to a velocity (B — > y/Jx^mrii B, with rrii the ion 
mass), P* is the total (magnetic plus kinetic) pressure and v is the viscosity. 
The Hall effect appears in the induction equation as an additional term pro- 
portional to the ion inertial length di which means that it is effective when the 
dynamical scale is small enough. In other words, for large scale phenomena 
this term is negligible and we recover the standard MHD equations. In the 
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opposite limit, e.g. for very fast time scales (<§; c^ 1 ), ions do not have time to 
follow electrons and they provide a static homogeneous background on which 
electrons move. Such a model where the dynamics is entirely governed by 
electrons is called the Electron MHD (EMHD) model (Kingsep et al., 1990; 
Shukla and Stenflo, 1999). It can be recovered from Hall MHD by taking 
the limits of small velocity V and large di. The Electron and Hall MHD ap- 
proximations are particularly relevant in the context of collisionless magnetic 
reconnection where the diffusion region develops multiscale structures corre- 
sponding to ion and electron characteristic lengths (Huba, 1995; Biskamp, 
1997). For example, it is often considered that whistler/EMHD turbulence 
may act as a detector for magnetic reconnection at the magnetopause (Cai 
et al., 2001). 



2.3 Three-dimensional inviscid invariants 

The three inviscid (u — rj — 0) quadratic invariants of incompressible Hall 
MHD are the total energy 

£ = ^/(V 2 + B 2 )dV, (6) 

the magnetic helicity 

H m = \j A-BdV, (7) 
and the generalized hybrid helicity 

H G = ^J(A + d,V) • (B + d t V x V) d V , (8) 

with A the vector potential (B = V x A). The third invariant generalizes the 
cross-helicity, H c — (1/2) /V • Be? V, which is not conserved anymore when 
the Hall term is present in the MHD equations. The generalized hybrid 
helicity can be seen as the product of two generalized quantities, a vector 
potential T = A + rfjV and a vorticity 17 = B + <ijV x V. The role played 
by the generalized vorticity is somewhat equivalent to the one played by the 
magnetic field in standard MHD (Woltjer, 1958). Indeed, both quantities 
obey the same Lagrangian equation, which is for fl, 

? = fi-VV. (9) 

dt 
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By applying the Helmholtz's law (see e.g. Davidson, 2001) to Hall MHD, we 
see that the generalized vorticity lines are frozen into the plasma (see e.g. 
Sahraoui et a!., 2003). The presence of the Hall term breaks such a property 
for the magnetic field which is however still frozen but only in the electron 
flow: the introduction of the electron velocity V e , with VxB-JxB/tie ~ 
V e x B, leads to 

f = B.VV., (10) 

which proves the statement. We will see below that the detailed conservation 
of invariants is the first test that the wave kinetic equations have to satisfy. 



2.4 Incompressible Hall MHD waves 

One of the main effects produced by the presence of the Hall term is that the 
linearly polarized Alfven waves, solutions of the standard MHD equations, 
become circularly polarized and dispersive (see e.g. Mahajan and Krishan, 
2005; Sahraoui et a!., 2005). Indeed, if we linearize equations around 
a strong uniform magnetic field B such that, 

B(x) = B e\\ +eb(x), (11) 

V(x) = ev(x), (12) 

with e a small parameter (0 < e C 1), x a three-dimensional displacement 
vector, and en a unit vector (|e||| = 1), then we obtain the following inviscid 
equations in Fourier space: 

k-v k = 0, (13) 

<9 f v k - zfc||B b k = e {-v • V v - VP* + b ■ V b} k , (14) 

<9 t b k — i/c||-BoUk — diB ki\ kxb k = e {— v- V b + b- V v — dt V x [(Vxb)xb]} k , 

(15) 

k-b k = 0, (16) 

where the wavevector k = ke^ = kj_ + fc||e|| (k = |k|, k± = |kj_|, |e^| = 1) and 
i 2 = —1. The index k denotes the Fourier transform, defined by the relation 

v(x) = / v(k)e*- x dk, (17) 
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where v(k) = = v^e luJt (the same notation is used for the magnetic 
field). The linear dispersion relation (e = 0) reads 



with 

We obtain the solutions 





{AdiB k\\k)uj- £gjfeg = 0, (18) 
Aie fc x(~ k ). (19) 

(20) 

where the value (±1) of s defines the directional wave polarity. In other 
words, we have skn > and is a positive frequency. The Alfven wave 
polarization A tells us if the wave is right (A = s) or left (A = — s) circularly 
polarized. In the first case, we are dealing with whistler waves, whereas 
in the latter case with ion cyclotron waves. We see that the transverse 
circularly polarized Alfven waves are dispersive and we note that we recover 
the two well-known limits, i.e. the pure whistler waves (u = skukdiBo) in 
the high frequency limit (kdi — > oo), and the standard Alfven waves (u = 
sk»B ) in the low frequency limit (kdi ~^ 0). The Alfven waves become 
linearly polarized only when the Hall term vanishes: when the Hall term is 
present, whatever its magnitude is, the Alfven waves are circularly polarized. 
Note that this situation is different from the compressible case for which the 
Alfven waves are elliptically polarized. As expected, it is possible to show 
(Hameiri et al., 2005; Sahraoui et al., 2005) that the ion cyclotron wave has 
a resonance at the frequency Udkn/k, where uo C i = Bo/di. Therefore, with 
such an approximation, only whistler waves survive at high frequency. 



2.5 Complex helicity decomposition 

Given the incompressibility constraints (fT3*J) and it is convenient to 

project the Hall MHD equations in the plane orthogonal to k. We will use 
the complex helicity decomposition technique which has been shown to be 
effective in providing a compact description of the dynamics of 3D incom- 
pressible fluids (Craya, 1958; Moffatt, 1970; Kraichnan, 1973; Cambon et al., 
1989; Lesieur, 1990; Waleffe, 1992; Turner,2000; Galtier, 2003; Galtier and 
Bhattacharjee, 2003). The complex helicity basis is also particularly useful 
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since it allows to diagonalize systems dealing with circularly polarized waves. 
We introduce the complex helicity decomposition 

h A (k) = h A = e e + *Ae $ , (21) 

where 

e e = e$ x e k , (22) 

e$ = rj ^7 , 23 

ley x e k \ 

and |eg(k)| = |e$(k)| = l. We note that (e^, h^, h^) form a complex basis with 
the following properties: 

h k A = h A k , (24) 

e fe xh A = -zAh A , (25) 

k-h A = 0, (26) 

h A -h A ' = 2 5_ A , A . (27) 

We project the Fourier transform of the original vectors v(x) and b(x) on 
the helicity basis: 

Vk =^w A (k)h A =5:w A h A , (28) 

A A 

b k =^s A (k)h A =5:s A h A . (29) 

A A 

We introduce expressions of the fields into the Hall MHD equations written 
in Fourier space and we multiply by vector h A k . First, we will focus on the 
linear dispersion relation (e = 0) which reads: 

d t Z s A = -iu s A Z s A , (30) 

with 

ZA=Z/ A +&B A , (31) 



Equation ()30|) shows that Z A are the "good" variables for our system. These 
eigenvectors combine the velocity and the magnetic field in a non trivial way 
by a factor £ A (with uj s k = — Bok\\£ A ). In the large scale limit (kdi — ► 0), 
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we see that £ A — > — s; we recover the Elsasser variables used in standard 
MHD. In the small scale limit (kdi — > oo), we have £ A — > —sdik, for A = s 
(whistler waves), or £ A — ► (— sdj/c) -1 , for A = — s. Therefore 2 A can be seen 
as a generalization of the Elsasser variables to Hall MHD. In the rest of the 
paper, we will use the relation 

2i = (a-*rKe-***. (33) 

where a% is the wave amplitude in the interaction representation for which 
we have, in the linear approximation, d t a s A = 0. In particular, that means 
that weak nonlinearities will modify only slowly in time the Hall MHD wave 
amplitudes. The coefficient in front of the wave amplitude is introduced in 
advance to simplify the algebra that we are going to develop. 



3 Helical wave turbulence formalism 
3.1 Fundamental equations 

We decompose the inviscid nonlinear Hall MHD equations (|14 )) -([15 )) on the 
complex helicity decomposition introduced in the previous section. Then we 
project the equations on vector h A k . We obtain: 

d t U A - iB k ]{ B A = (34) 
/ E (*V4, -B Ap B Aq )(k-h^)(h^ -hk A ) 5 m , k dpdq, 



and 



le 
~2 



d t B A - iB Q k\pi A + iAdiB Q k\\kB A = (35) 

-4 I E( w A^A 9 - w A^A P )( k - h ? p )K q - h k A )W^p^q 

A P ,A ? 

+ 

[kAd t J2 i3 Ap S A J(q-h k A )(h^-h^)-(q-h^)(h^-h k A )]Wrfprfq, 

J A A 



where 6 pq ^ = S(p + q — k). The delta distributions come from the Fourier 
transforms of the nonlinear terms. We introduce the generalized Elsasser 
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variables a s A in the interaction representation and we find: 




AApAq 



dta s A = - 



1 W dp dq , 



(36) 



where 




i-q 2 
a - a- 



7(q-^)(hJ 



<*) + 



(37) 



k p q 



SASA g — SA P SA 9 



(qh^)(h^.h£), 



£1 — £a 



and 



(38) 



Equation (pffij) is the wave amplitude equation from which it is possible to 
extract some information. As expected we see that the nonlinear terms are of 
order e. This means that weak nonlinearities will modify only slowly in time 
the Hall MHD wave amplitude. They contain an exponentially oscillating 
term which is essential for the asymptotic closure. Indeed, wave turbulence 
deals with variations of spectral densities at very large time, i.e. for a non- 
linear transfer time much greater than the wave period. As a consequence, 
most of the nonlinear terms are destroyed by phase mixing and only a few of 
them, the resonance terms, survive (see e.g. Newell et al., 2001). The expres- 
sion obtained for the fundamental equation (J36|) is usual in wave turbulence. 
The main difference between problems is localized in the matrix L which is 
interpreted as a complex geometrical coefficient. We will see below that the 
local decomposition allows to get a polar form for such a coefficient which 
is much easier to manipulate. From equation (jSHJ) we see eventually that, 
contrary to incompressible MHD, there is no exact solutions to the nonlinear 
problem in incompressible Hall MHD. The origin of such a difference is that 
in MHD the nonlinear term involves Alfven waves traveling only in opposite 
directions whereas in Hall MHD this constrain does not exist (we have a 
summation over A and s). In other words, if one type of wave is not present 
in incompressible MHD then the nonlinear term cancels whereas in incom- 
pressible Hall MHD it is not the case. The conclusion reached by Mahajan 
and Krishan (2005) is therefore not correct: the condition U\ = — £a<Sa does 
not correspond to a nonlinear solution of incompressible Hall MHD. 
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Figure 1: Sketch of the local decomposition for a given wavevector p. 
3.2 Local decomposition 

In order to evaluate the scalar products of complex helical vectors found in 
the geometrical coefficient (|37|h it is convenient to introduce a vector basis 
local to each particular triad (Waleffe, 1992; Turner, 2000; Galtier, 2003; 
Galtier and Bhattacharjee, 2003). For example, for a given vector p, we 
define the orthonormal basis vectors, 

6«(p) =n, 

6( 2 )(p) =e p xn, (39) 
6( 3 )(p) =e p , 

where e p = p/|p| and 
n 

We see that the vector n is normal to any vector of the triad (k,p,q) and 
changes sign if p and q are interchanged, i.e. ri(k, q ,p) = — n(k, P , q )- Note that 
n does not change by cyclic permutation, i.e. , n(k, q , P ) = n( q , p ,k) = n( p k, q ). 
A sketch of the local decomposition is given in Fig. ^ We now introduce the 
vectors 

~Mp) = 2J* = 6( 1 )(p) + zA p 6( 2 )(p) , (41) 
and define the rotation angle $ p , so that 

cos$ p =n-e e (p), 
sin $ p = n ■ e$(p) 

14 



p x k 
|p x k| 



q x p 

lq x p| 



k x q 
Ik x q| 



(40) 



(42) 



The decomposition of the helicity vector hp p in the local basis gives 

h£p = e iAp ^ . (43) 
After some algebra we obtain the following polar form for the matrix L, 



L ss p s i = (44) 
k p q 

SA SA K 

+ 

a^c-eo^CAA.+cos^)]. 

The angle ipk refers to the angle opposite to k in the triangle defined by 
k = p + q (sin^ = n • (q x p)/|(q x p)|). To obtain equation (|44|) . we have 
also used the well-known triangle relations 

sinip k sinV'p sin^ 

— — = = . (45) 

k p q 

Further modifications have to be made before applying the spectral formal- 
ism. In particular, the fundamental equation has to be invariant under inter- 
change of p and q. To do so, we introduce the symmetrized and renormalized 
matrix M: 



M° s '> s i =idi — -p- q — ^ + ) . (46) 



kpq ' £a 2 6a 9 ~€a p \ k pq ' k qp 
Finally, by using the identities given in Appendix [XJ we obtain: 



C"9 _ £*P AApAg 
*± CLiJ a a _ C a C a K X) q 



A p ,A q 



where 



M SSpSq = (48) 
k p q 



e *(A* fc+ WA,*,) A A, A, ^ (Afc + A„ P + A^) (1 - £a ' V' V' 2 ) 
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The matrix M possesses the following properties (* denotes the complex 
conjugate), 



AApA, 

k p q 



) 



M 



A-Ap-A 9 



k p q 



M 



—k — p — q : 



A A p A q 



(49) 



AApAq 

M ss p s q 

k p q 



AAgAp 

M ss i s p 

k q p ' 



(50) 



AApAg 

M ss p s q 

k p q 



AqApA 
M s q s P s 

q p k 



(51) 



AApAq 

k p q 



ApAAq 
M SpSSq 

p k q 



(52) 



Equation (}4Tj) is the fundamental equations that describe the slow evolution 
of the Alfven wave amplitudes due to the nonlinear terms of the incom- 
pressible Hall MHD equations. It is the starting point for deriving the wave 
kinetic equations. The local decomposition used here allows us to represent 
concisely complex information in an exponential function. As we will see, it 
will simplify significantly the derivation of the wave kinetic equations. 

From equation (J47)) we note that the nonlinear coupling between helicity 
states associated with wavevectors, p and q, vanishes when the wavevectors 
are collinear (since then, sin-^fc = 0). This property is similar to the one 
found in the limit of EMHD (Galtier and Bhattacharjee, 2003). It seems 
to be a general property for helicity waves (Kraichnan,1973; Waleffe, 1993; 
Turner, 2000; Galtier, 2003). Additionally, we note that the nonlinear cou- 
pling between helicity states vanishes whenever the wavenumbers p and q 
are equal if their associated wave and directional polarities, A p , A q , and s p , 
s q respectively, are also equal. In the case of whistler (EMHD) waves, for 
which we have A = s (right circularly polarized), this property was already 
observed (Galtier and Bhattacharjee, 2003). Here we generalize this finding 
to right and left circularly polarized waves. In the large scale limit, i.e. when 
we tend to pure incompressible MHD, this property tends to disappear. For 
pure incompressible MHD, where Alfven waves are linearly polarized, this 
is not observed anymore. As noticed before, the nature of the polarization 
seems to be fundamental. 

We are interested by the long-time behavior of the Alfven wave ampli- 
tudes. From the fundamental equation (J47|) . we see that the nonlinear wave 
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coupling will come from resonant terms such that, 



k = p + q, 

(53) 



The resonance condition may also be written: 

SA SA P ?A 9 SA ?A 9 ?A P 

% h h 



(54) 



As we will see below, the relations (jHlj) are useful in simplifying the wave 
kinetic equations and demonstrating the conservation of ideal invariants. In 
particular, we note that we recover the resonance conditions for whistler 
waves by taking the appropriate limit. 



3.3 Dynamics and wave kinetic equations 

Fully developed wave turbulence is a state of a system composed of many 
simultaneously excited and interacting nonlinear waves where the energy dis- 
tribution, far from thermodynamic equilibrium, is characterized by a wide 
power law spectrum. This range of wavenumbers, the inertial range, is gen- 
erally localized between large scales at which energy is injected in the system 
and small dissipative scales. The origin of wave turbulence dates back to the 
early sixties and since then many papers have been devoted to the subject 
(see e.g. Hasselmnan, 1962; Benney and Saffman, 1966; Zakharov, 1967; Ben- 
ney and Newell, 1969; Sagdeev and Galeev, 1969; Kuznetsov, 1972; Zakharov 
et al.; 1992; Newell et al., 2001). The essence of weak wave turbulence is the 
statistical study of large ensembles of weakly interacting dispersive waves 
via a systematic asymptotic expansion in powers of small nonlinearity. This 
technique leads finally to the derivation of wave kinetic equations for quan- 
tities like the energy and more generally for the quadratic invariants of the 
system studied. Here, we will follow the standard Eulerian formalism of wave 
turbulence (see e.g. Benney and Newell, 1969). 

We define the density tensor g A (k) for an homogeneous turbulence, such 
that: 

K(k) a A ',(k')> = q s A (k) 5(k + k') 5 AA , 5 SS , , (55) 

for which we shall write a "closure" equation. The presence of the delta 5aa' 
and 5 ss i means that correlations with opposite wave or directional polarities 



17 



have no long-time influence in the wave turbulence regime; the third delta 
distribution 5(k + k') is the consequence of the homogeneity assumption. 
It is strongly linked to the form of the frequency uj a (see the discussion in 
Section Details of the derivation of the wave kinetic equations are given 
in Appendix El We obtain the following result: 



%l(k) 



vre 2 


f E 

A p ,Ag 
s p-, s q 


AdjBlJ 


( < 









sin ipj. 



(Ak+A p p+A q q) 2 ( 1 - 



(56) 



C(p) 



a/ 2 / 



?A( k ) ?A p (p) ?£(<j) ^(^fe,p?) $k,pq dp dq . 

Equation (JHH)) is the main result of the helical wave turbulence formalism. It 
describes the statistical properties of incompressible Hall MHD wave turbu- 
lence at the lowest order, i.e. for three-wave interactions. 



4 General properties of Hall MHD wave tur- 
bulence 

4.1 Triadic conservation of ideal invariants 

In section 12 .3\ we have introduced the 3D ideal invariants of incompressible 
Hall MHD. The first test that the wave kinetic equations have to satisfy is 
the detailed conservation of these invariants, that is to say, the conservation 
of invariants for each triad (k, p, q). Starting from definitions (JO}-©? we 
note that the Fourier spectra of the ideal invariants are: 

£(k)=E(l+£r 2 )<?A(k), (57) 

A,s 

4(k)^T#). (58) 
# G (k) = E^T-9A(k). (59) 

A,s K 
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We will first check the energy conservation. From expression (JSBJ), we 



find: 



7r e 



AdfBQ J a ,a7aj V ^ 



d t E(t) =d t [ E(k) dk = d t [j2 4(k) rfk = 

A,s 



(60) 

I fell 



^A 



e A (k) e £(p) e£(q) 



9a ( k ) 9a!, (p) 9a' (o) h m dk dp dq . 



Equation (|50jl is invariant under cyclic permutations of wavevectors. That 
leads to: 

d t E(t) = (61) 



7T e 



12 d? R 2 



E 



/ sin ifck 



i ^0 J A, A 



A \ ^ 



(Ak+A pP +A q qY [ i - er^re 



I fell 



3,Sp,Sq 



P'l 



^a; 



"A 



e A (k) c *(p) e £(q) 



9l( k ) 9a p (p) ?a* (q) 5 (^fc, P g) rfkrf p rf q- 



Total energy is conserved exactly on the resonant manifold since then Qk,pq = 
0: we have triadic conservation of total energy. 

For the other invariants, it is straightforward to show that the difference 
between them is conserved by the wave kinetic equations. With relation 
), we obtain: 



d t (H m (t) - Hq (t) ) = dtj (H m (k) - H G (k)) dk 



n e 



E 



'sin^ 



(Ak+A pP +A qq y i - er c £ 



(62) 



fen 



^A 



9a ( k ) 9a* (p) ?A ? (q) $(ttk, P q) S k , Pq dkdp dq. 
Equation (|5Uj) is also invariant under cyclic permutations of wavevectors. 



9A;iqJ 



Then one is led to 



d t (H m (t) - H G {t)) 



(63) 
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i + erv aw U + c 2 / ^ (p) vi + C / #(<o 

9A ( k ) 9a* (p) ?A ? (q) <K^fc,w) rfkrfp dq, 

which is exactly equal to zero on the resonant manifold: we have triadic 
conservation of the difference between magnetic and global helicity. The 
last sufficient step would be to show detailed conservation for one of the 
two helicities. Unfortunately it is not so trivial and we will not give the 
proof here. However, we note already that magnetic helicity is conserved for 
equilateral triangles (k = p = q), i.e. for strongly local interactions. 



(h-p\\-n) 



4.2 General properties 

From the wave kinetic equations (JEEJ), we find several general properties. 
Some of them can be obtained directly from the wave amplitude equation 
f!47|) as explained in Section ET2l First, we observe that there is no coupling 
between Hall MHD waves associated with wavevectors, p and q, when the 
wavevectors are collinear (sin-?/^ = 0). Second, we note that there is no 
coupling between helical waves associate with these vectors whenever the 
magnitudes, p and q, are equal if their associated polarities, s p and s q in 
one hand and, A p and A q on the other hand, are also equal (since then, 
£a. 9 — £lp = 0)- These properties hold for the three inviscid invariants and 
generalize what was found previously in EMHD (Galtier and Bhattacharjee, 
2003) where we only have right circularly polarized waves (A = s). It seems to 
be a generic property of helical wave interactions (Kraichnan,1973; Waleffe, 
1993; Turner, 2000; Galtier, 2003). As noted before, this property tends to 
disappear when the large scale limit is taken, i.e. when we tend to standard 
MHD. Third, it follows from the previous observations that a strong helical 
perturbation localized initially in a narrow band of wavenumbers will lead 
to a weak transfer of energy, magnetic and global helicities. Note that these 
properties can be inferred from the fundamental equation ({17}) as well. 
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4.3 High frequency limit of electron MHD 

In the present section we shall demonstrate that the small scale limit (dik — ► 
+00) of the wave kinetic equations (JSTIj) tends to the expected equations of 
electron MHD when only right (A = s) circularly polarized waves are taken 
into account. We recall that, in the small scale limit, £| — > —sdik and 
^ s — > s/dik. One obtains 

d,qt(k) = (64) 

^ / E (nr ) 2 («* + w + w>" ( M ^ M ) 2 

_«ArAc|j g£(p) g£(q) - s p pp|| g s s (k) g»«(q) - s q qq\\ q s s {k) q%(p)] 5(Q k , pq ) 8 km dpdq : 

where flk, pq = B di(sknk — s p p\\p — s q q\\q). The kinetic equations found have 
exactly the same form as in Galtier and Bhattacharjee (2003) (where by 
definition uj s s = su^, B = 1 and dj = 1) who derived a wave turbulence 
theory for electron MHD (see the discussion in Section IS)) . In particular, this 
means that by recovering the electron MHD description from the present 
Hall MHD theory, we recover all the properties already found for whistler 
wave turbulence (anisotropy, scaling laws, direct cascade...). We will come 
back to that point in Section |S| Note finally that there is a strong analogy 
between such a limit and wave turbulence in rapidly rotating flows (Galtier, 
2003; Bellet et al., 2005; Morize et al., 2005). The physical reason is that 
in both problems (i) there is a privileged direction, played by the rotating 
axis or the magnetic field B , (ii) there are dispersive helical waves, called 
inertial waves for rotating flows, and (iii) the wave frequencies are not very 
different, i.e. the inertial wave frequency is proportional to k\\/k. In the past, 
comparisons have been made between incompressible rotating turbulence and 
magnetized plasmas described by incompressible MHD. The results obtained 
in the framework of wave turbulence show clearly that the comparison is 
much more relevant if one considers electron MHD plasmas. Indeed, a small 
nonlinear transfer is found along the privileged direction for rotating and 
electron MHD turbulence whereas this kind of transfer is strictly forbidden 
in MHD turbulence. 

4.4 High frequency limit of ion MHD 

We may also be interested by the small scale limit (dik — > +00) of the wave 
kinetic equations (J56)) when only left (A = — s) circularly polarized waves 
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are taken into account. We decide to call this limit the ion MHD (IMHD) 
approximation following the example of the electron MHD. We have, for such 
a limit, £z| — ► sdik and £f_ s — ► —s/dik. One finds 



E 



/ smip k 

V k 



(sk + s p p + s q q) 



2 / s qQ S pP 



pqk\\ 



(65) 



sk\\kp q 5(£l k , P q) S k , pq 



s P P\\ 



q s _ s (Vq- q s ,(q) 



M«^(k)^(p) 



with Qk, P q = {Bq/ di){sk\\jk — Sppn/p — s q q»/q). No existing theory has been 
developed for such a limit. We will see in Section |S] that it is possible to 
extract from the master equations of ion cyclotron turbulence some exact 
properties. 



4.5 Low frequency limit of standard MHD 

This section is devoted to the low frequency limit, i.e. the large scale limit, of 
the wave kinetic equations (|56p of Hall MHD. The MHD limit is somewhat 
singular for Hall MHD. Indeed the large scale limit does not tend to the 
expected wave kinetic equations that were derived first by Galtier et al. 
(2000). The subtle point resides in the kinematics: for pseudo-dispersive 
MHD waves, the definition of the density tensor q%(k) that is used 
for Hall MHD is not valid anymore. The reason is that in the large scale 
limit the polarization A does not appear anymore in the frequency, which is 
uj = sk\\B Q . Then the kinematics tells us that the definition for the density 
tensor is 

K(k) 4(k')> = q s A i *(k + k') 5 SS , , (66) 

where the condition A = A' is not necessary satisfied. This means that the 
large scale limit of equations (J56|) leads to MHD wave kinetic equations in 
the particular case where helicity terms are supposed to be absent and where 
equality between shear-Alfven and pseudo-Alfven waves energy is assumed. 
Indeed, helicity terms involve quantities for which A = —A' and the total 
energy involves only terms for which A = A'. However, energies for shear- 
Alfven waves and pseudo-Alfven waves involve terms with different polarities. 
For shear-Alfven waves, we have: 

v k - sb k = Z;h+ + Zlh k = (Z' + - Z s _) % e„ , (67) 
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and for pseudo-Alfven waves: 



v k - sb k = Z|h+ + Z s _h k = {Z s + + Zt) e e . (68) 
Energies associated to shear- and pseudo-Alfven waves are respectively: 

(z s + - z s _)(z s + - z s _y = \z s + \ 2 + \zt\ 2 - z s + (z s _y - z s _(z s + y , (69) 

and 

(z» + z s _){z s + + zty = \z%\ 2 + \z s _\ 2 + z s + (z s _y + z s _(z s + y . (70) 

We see clearly that if, in the wave kinetic equations, we only take into ac- 
count (quadratic) terms with the same polarizations (A) then it is equivalent 
to assume equality between shear- and pseudo-Alfven wave energies. We 
will see below that the large scale limit of the wave kinetic equations (|56|) 
of Hall MHD tends to the expected MHD counterpart when the previous 
assumptions about helicities and energies are satisfied. A derivation of the 
wave kinetic equations is given in Appendix The result is given for the 
density tensor q^. Contrary to Hall MHD, and actually to any problem 
in wave turbulence, principal value terms appear for incompressible MHD. 
The reason of the presence of principal value terms is linked to the nature of 
Alfven waves which are pseudo-dispersive. 

Further comparisons between the results in Appendix and the wave 
kinetic equations obtained by Galtier et al. (2000) are of course possible (see 
e.g. the discussion in Section [HJ) but, for simplicity, we prefer to focus our 
attention to the case where the density tensor is symmetric in A. Therefore 
we start our analysis with the general kinetic equation ()56|) and take the 
large scale limit (dik —>■ 0) for which we have, at the leading order, £^ ~~ * 
—s — Adik/2. After some simplifications, we arrive at: 

dtql(k) = (71) 

i£ / E (nr ) 2 ( Ak + A P p + A ^) 2 ( skA + v?a p + s qq k q f fej^j 2 

sp,sq 

sk \\ ( sfc ||9A p (p)9A,(q) - s pWA( k )9A,(q) - s ?9||?A( k )9Ap(p)) 8(n k , pq )5 kiPq dpdq, 

where flk, pq = B (sk\\ — s p p» — s q qn). Note that we only have a nonlinear 
contribution when the wave polarities s p and s q are different. We recover 
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here a well-known property of incompressible MHD: in such a limit, we only 
have nonlinear interactions between Alfven waves propagating in different 
directions. One expands the summation over the directional polarities s p 



and s q , and finds 




«(k) = (72) 
(Ak + A p p + A q q) 2 (Ak + A p p - A q q) 2 



Qa p (p) -<?l( k )J S(q\\)5 k<pq dpdq. 

This result is exactly the same as in Appendix \D\ when MHD wave kinetic 
equations are considered in the particular case where only terms symmetric 
in A are retained, i.e. terms like q AA . Therefore under these assumptions the 
MHD description does not appear like a singular limit for Hall MHD. (Note 
that this continuity in the description is a priori not so trivial if we remember 
that as soon as the Hall term is included in the standard MHD equations, 
whatever its magnitude is, it changes the polarity of the Alfven waves which 
become circularly polarized.) 

The comparison with the wave kinetic equations derived by Galtier et 
al. (2000) (see equations (26)) is not direct since here the problem has been 
decomposed at the origin on a complex helicity basis. The main signature of 
this decomposition is the dependency on the wave polarity A. In spite of this 
difficulty, a common point is clearly seen in the presence of the delta S(q\\) 
which arises because of the three-wave frequency resonance condition. This 
means that in any triadic resonant interaction, there is always one wave that 
corresponds to a purely 2D motion (gii = 0) whereas the two others have 
equal parallel components {jpu = kn). The direct consequence is the absence 
of nonlinear transfer along B , a result predicted earlier by several authors 
(see e.g. Montgomery and Turner, 1981; Shebalin et al., 1983). In other 
words, we have a two-dimensionalization of the Alfven wave turbulence (see 
also e.g. Ng and Bhattacharjee, 1997; Lithwick and Goldreich, 2003). 



5 Master equations of Hall MHD turbulence 

5.1 General case 

In order to extract further information about Hall MHD turbulence, we are 
going to write the expression of the spectral density g A (k) in terms of the 
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invariants. In practice, we need to add a fourth variable, Ed, which has been 
chosen to be the difference between the kinetic and magnetic energy, namely 



^(k) = E(^ s - 1) <ZA(k) ■ (73) 

Note that contrary to the MHD case, wave turbulence in Hall MHD allows 
to have a departure from equipartition between kinetic and magnetic energy 
and therefore a non trivial value for Ed- We will see, in Section that 
this property is fundamental to get non trivial nonlocal interactions between 
waves. It is possible to inverse the system q\(E, Ed, Hm, Hq); one finds 

£00 = (74) 

(a 2 + a s2 )[(a 2 - g) - (a 2 + mm + 2Ak(g i H m (k) - Ham 

4 (a 4 - cr 4 ) 

The introduction of expression (|7ijl into (|56|) leads to the wave kinetic equa- 
tions for E, Ed, Hm and Hq- We will not make such a lengthy development 
and we will rather focus on the master equations which drive the Hall MHD 
turbulence. Indeed, because of the presence of the factor A in expression 
(174)1 . we see that the nonlinear terms with different polarities will not play 
the same role. In particular, for the wave kinetic equations of energies, only 
the interactions between nonlinear terms involving either energies or helici- 
ties will give a contribution. For the wave kinetic equations of helicities, we 
only have contributions from nonlinear terms involving energies and helici- 
ties. Thus the energy equations are the master equations driving Hall MHD 
turbulence. In other words, this means that an initial state with zero helic- 
ity will not generate any helicity at any scale. However, an initial state of 
zero helicity does not preclude the development of energy spectra. Assum- 
ing H m (k) = and i?c(k) = 0, the master kinetic equations of Hall MHD 
turbulence, written for the kinetic energy E v and the magnetic energy E , 
are 
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fcn 



, 2 




In Appendix [Ejthe equivalent kinetic equations for the variables E and Ed are 
given. Below we describe the small scale (whistler and ion cyclotron waves) 
and large scale (pure Alfven waves) limits of such a description. In particular, 
we will see that the role played by the kinetic and magnetic energies may be 
very different and that the master equations can be simplified further. 



5.2 Master equations in the limit of Whistler wave tur- 
bulence 

In the small scale limit (dik — > +00), one can distinguish between the whistler 
branch and the ion cyclotron branch. For right circularly polarized wave 
(A = s) we have — > l/d 2 k 2 , whereas for left circularly polarized wave 
(A = — s), — ► d 2 k 2 . Then in the limit of whistler wave turbulence 
(EMHD limit), we obtain 

with flk, Pq = Bq di(sk\\k—s p p\\p—s q q\\q). We see clearly that the kinetic energy 
is not a relevant quantity in such a limit since the factor l/dfk 2 damps the 
nonlinear contributions. This is not in contradiction with what we know 
if we remember that the velocity in the original Hall MHD equations is a 
combination of the ion and electron velocities. In the EMHD limit, i.e. in 
the small spatio-temporal scale limit, this velocity tends to be small: ions 
do not have time to follow electrons and they provide a static homogeneous 
background on which electrons move. Since the magnetic field is frozen 
into the ideal electron flow (see Section 12.31 for the proof) we understand 




sfc||fcs p P|| P E"(q) E B (p)-E B (k) 
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why the magnetic energy is the dominant contribution to the wave kinetic 
equations. The wave kinetic equation found for the magnetic energy have, 
at leading order, exactly the same form as in Galtier and Bhattacharjee 
(2003) (where u% = su>k and B = 1) but they are more general in the 
sense that we have at our disposal not only an equation for the magnetic 
field, the dominant contribution, but also for the velocity that behaves like 
the magnetic field. The consequence is that we recover all the properties 
already known. In particular, in such a regime the energy spectrum follows 
the power law E(k±,k\\) ~ k ± 5 ^ 2 k^ [k\\ is assumed positive) which is an 
exact solution of the wave kinetic equations. 



5.3 Master equations in the limit of ion cyclotron wave 
turbulence 

For the ion cyclotron branch (A = — s), the small scale limit (dik — > +oo) 
gives 



7r e 
8£ 



E 



dl EV{k) \~ 

'^)\sk +SpP+Sqq r(^ 



(77) 



dfk 1 sk\\s p p\\pq 



k 



E v (q) E v (p)-E v (k) 5(n ktPq )5 ktPq dpdq, 



with Qk, pq = (Bo/di)(sk\\/k — s p p\\/p — s q q\\/q). In such a limit, we are able 
to describe both the kinetic and the magnetic energies but we see that the 
dominant quantity is the kinetic energy since the factor d\k 2 amplifies the 
nonlinear contributions in the corresponding equation (see e.g. Mahajan and 
Krishan, 2005). The physical reason is that in such a limit, we are mainly 
dealing with ions whereas the magnetic field is frozen into the ideal electron 
flow (see Section 12. 3|) . Like for the EMHD limit, it is possible to derive 
the exact power law solutions. First of all, we note that for strongly local 
interactions the resonance condition writes 



Sq ~ S 



P\\ 



S — S r 



n 



(78) 



which leads, as for the whistler waves case, to anisotropic turbulence. The 
same analysis for strongly nonlocal interactions leads to the same conclusion. 
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Thus, we may write the wave kinetic equations in the limit k± 3> fen. For an 
axisymmetric turbulence E v (\t±, k\\) = E v (k±, k\\)/2'jrk±; we find 



d t E v (k ± ,k l{ 



(79) 



Y^J sin Vv ( sk ± + s pP± + s q1±) 2 [ 



s P p± - s q q± 



sk 



\\ s pP\\ 



E V (q±,q\\) [k ± E v (p ± , Pll ) -p^E v (k Ll k\\ 



k\\ P± q± 
8(ttk, pq ) f>k p p m dp ± dq ± dp\\dq\i 



where flk, pq = (B /di)(sk\\/k± — SpP\\/p± — s q q\\/q±). Note that in such a limit 
the well-known triangle relations are: smip k± /k± = sin ip p± /p±_ = sin ip q± /q±. 
The wave kinetic equation ()79)1 is now symmetric enough to apply a con- 
formal transformation, called the Kuznetsov-Zakharov transformation. This 
transformation, a two-dimensional generalization of the Zakharov transfor- 
mation, has been applied in various anisotropic problems (Kuznetsov, 1972; 
Balk et al. , 1990; Kuznetsov, 2001; Galtier and Bhattacharjee, 2003; Galtier, 
2003). The bihomogeneity of the collision integrals in the wavenumbers k± 
and k\\ allows us to use the transformation 



P± 
Q± 
P\\ 

n 



k±/p± , 

k±q±/p± 
k f\/P\\ > 

hn/pw ■ 



(80) 



We search for stationary solutions in the power law form E(k±, k\\) ~ kj_ n k^ m . 
(We will only consider positive parallel wavenumber.) The new form of the 
collision integral of equation (j79jl . resulting from the summation of the inte- 
grand in its primary form and after the Kuznetsov-Zakharov transformation, 
is 

d t E v {k ± ,k {] ) = (81) 

/ \ 2 

^ torn. on, \ 

sxnip q± (sk x +s p p±+s q q ± 
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)2 ( SpPL ~ s q q± 
\ k\\ p ± q ± 



pq ) 



p± 

k i 



— n— 1 




P±\ 

kj 



2n-5 




2m- 1 N 



dp±_ dq± dp\\ dq\\ . 
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The above collision integral vanishes for specific values of n and m. The exact 
power law solutions correspond to these values. There are two different kind 
of solutions. The fluxless solution, also called the thermodynamic equilibrium 
solution, correspond to the equipartition state for which the flux of energy is 
zero. For this case, we have 

n = — 1 , 



This result can easily be checked by direct substitution in the original wave 
kinetic equation. The most interesting solution of the wave kinetic equa- 
tion (J81|) is the one for which the flux is non-zero and finite. The exact 
solution is called the Kuznetsov-Zakharov-Kolmogorov (KZK) spectrum and 
corresponds to the values, 

n =5/2, 

i in (83) 
m = 1/2. v ' 

In other words, the KZK solution scales as 

E v {k ± ,k\{) ~ £^ 5/ \ 1/2 . (84) 

We note that the power law scaling found is the same as the one for whistler 
wave turbulence. However here we are dealing with the kinetic energy not 
the magnetic one. A necessary condition for the realizability of the KZK 
spectra is that the turbulence be local in the sense that the behavior of the 
turbulence is determined primarily (but not only) by interaction between 
wave packets of comparable spatial scale (see e.g. [I]). To check a posteriori 
the validity of the solutions, we need to determine the domain of locality, 
i.e., the domain where the collision integral converges. In practice, we check 
that the contribution of nonlocal interactions do not lead to a divergence of 
the collision integral. Here, the condition of locality is automatically satisfied 
since the anisotropic limit introduces in the problem a cut-off which prevents 
any infra-red divergence of the collision integral. This characteristic is also 
observed for internal gravity waves (Caillol and Zeitlin, 2000), inertial waves 
(Galtier, 2003) and whistler waves (Galtier and Bhattacharjee, 2003). The 
other consequence of the existence of such a cut-off is that it is not possible to 
evaluate precisely the Kolmogorov constant, i.e. the prefactor of the spectra. 
However, the sign of the energy transfer can be computed for a reasonable 
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range of cut-offs. We observe a positive sign which means that we have a 
direct energy cascade. Note that this information cannot be obtained through 
a heuristic reasoning but only with a rigorous analysis. 



5.4 Master equations in the limit of pure Alfven wave 
turbulence 

In the large scale limit (dik — > 0) for which terms like tend to 1, we 
note that an equipartition between the kinetic and magnetic energies may 
be obtained since their wave kinetic equations tend to be identical: if ini- 
tially there is equipartition (E v (k) = E B {k) = E{k)/2) then at any time 
the equipartition will be conserved. The large scale limit corresponds to the 
standard MHD approximation for which the equipartition is in fact automat- 
ically satisfied by the kinematics. The large scale limit of expression (|75|) . 
for a state of equipartition, gives 



d t E(k) = (85) 



2 / • / \ 2 

7r e f ^ / sin ip: ^ 
256 



/ E ( ) ( Ak + A pP + V) 2 {sAk + s p A pP + s q A q qf 



, s p, s q 



2 

| sk\\s p p\\ E(q) [E(p) - E(k)} 5{tt km ) 5 km dp dq . 



with Vt kpq = B (sk\\ — s p p\\ — s q q\\). As expected, only nonlinear interactions 
between Alfven waves with different directional polarities s p = — s q will con- 
tribute to the dynamics. After such a consideration and an expansion over 
the polarities, one finds 

d t E(k) = (86) 

2 

J sin 2 ip k (1 + cos 2 ip g ) p 2 E(q) [E(p) - E(k)} S(q { \) 5 k>pq dp dq , 
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where we have used the identities 



£ {Ak + A p p + A q q) 2 {Ak + A pP - A q qf = (87) 

A,A p ,Aq 

{{k + p) 2 - q 2 ) 2 + {{k - p) 2 - q 2 ) 2 ] = 32k 2 p 2 {l + cos 2 ij> q ) , 
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and the well known triangle relations. The angle ijj q refers to the angle 
opposite to q in the triangle defined by k = p + q. (The same result can 
actually be found directly from equation ([72)1: the inversion of the system 
made in section I5~T1 is thus compatible with such a limit.) As noted before, 
the large scale limit of Hall MHD leads to the particular case where helicities 
are absent and where shear- and pseudo-Alfven wave energies are the same. 
Equation (J86j) can be recovered from the general kinetic equations obtained 
by Galtier et al. (2000) when the same assumptions are made. The proof 
is given in Appendix [FJ Therefore, all the properties previously derived are 
recovered (scaling laws, direction of the cascade, Kolmogorov constant...). In 
particular, the energy spectrum follows the power law E(k±, k») ~ kj 2 f(k\\), 
where / is an arbitrary function that is due to the dynamical decoupling of 
parallel planes in Fourier space. This is an exact solution of the wave kinetic 
equations. 



6 Anisotropic heuristic description for Hall 
MHD 

The study of wave turbulence in Hall MHD is a difficult task. We have seen 
above that for three specific limits we are able to find the exact power law 
energy spectra. Two of them (Alfven and whistler wave turbulence) were 
already known and a heuristic description was given by Galtier et al. (2000) 
and Galtier and Bhattacharjee (2003). We shall derive here a generalized 
heuristic description able to recover the essential physics underlying the KZK 
spectra including the ion cyclotron wave turbulence. However, note that 
several phenomenologies may be used to take into account the various regimes 
involving a non trivial balance between propagation and nonlinear effects like 
for standard MHD (see e.g. Zhou et al., 2004). 

We introduce anisotropy in our description by considering that 

k « k± > k\\ . (88) 

The primary variables used in the formalism are the generalized 
variables Z s k . Thus the nonlinear time built on the generalized 
variables is 

1 
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Note here that Z S K has a dimension of a velocity. In other words, it is not 
taken in Fourier space as it was introduced in Section 1231 In a similar way, 
we find the following Hall MHD wave period 

T ~~^r^ka- (90) 

We introduce now the mean rate of energy dissipation per unit mass II. Con- 
trary to the electron MHD case, we do not have to renormalize this quantity 
since it is automatically taken into account by the generalized Elsasser vari- 
ables. Then we have 

E E(k±,k\\)k±k\\ . s 

n v \y_LA ^ ( 91 ) 

Ttr T tr 

where the transfer time Tt r has the usual form given by the wave kinetic 
equation 

Ttr ~ 7iV L — ; (92) 
Tw 

it gives 

n~ 1 'J'V ■ (93) 

To proceed further, we note that 

(zi) 2 ~(a-a A s ) 2 < 2 , (94) 

and 

£~(l + a s2 )4 2 - (95) 

In relations (JHlj) and (jHSJ), both and a s A are written in the physical space, 
not in Fourier space. By using the relationships given in Appendix we 
finally obtain 



E{k ± , k\\) ~ ^n^o kl 2 k~ 1/2 {l + kid 2 )- 1 / 4 . (96) 

The heuristic prediction proposed here is able to describe anisotropic tur- 
bulence for the three different limits discussed above. We recover, in the 
small scale limit (k±di — > oo), the expected scaling law for whistler as well 
as ion cyclotron wave turbulence and, in the large scale limit (k±di — > 0), 
the Alfven wave turbulence scaling law (since then the parallel wavenumber 
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is a mute variable). The prediction is given for the total energy, therefore 
for the previous small scale limits one needs to consider only the magnetic or 
the kinetic energy respectively. One may understand this point by looking at 
equation (j3*Tj) where the relation between the generalized Elsasser variables, 
the velocity and the magnetic field is given. In the small scale limit, we have 
for whistler waves — > —sdik and for ion cyclotron waves £j[ — > {—sdik)^ 1 
(whereas in the large scale limit £^ — > — s); thus we see that either the ki- 
netic or the magnetic field will dominate. Note that we ignore anisotropy 
and assume k± ~ ku ~ k, we arrive at the scaling 

E(k) ~ yTLB^A;- 3 / 2 (l + k 2 dl)~ l,A (97) 

for the one- dimensional isotropic spectrum from which one recognizes the 
Iroshnikov (1963) and Kraichnan (1965) prediction for MHD. 

To summarize our finding we propose the picture given in Fig. [21 It 
is a sketch of the power law energy spectrum predictions for perpendicular 
wavenumbers at a given fen. The energy spectrum of Hall MHD is character- 
ized by two inertial ranges - the exact power law solutions of the wave kinetic 
equations - separated by a knee. The position of the knee corresponds to 
the scale where the Hall term becomes sub/dominant, i.e. when k±di ~ 1. 
The heuristic prediction tries to make the link continuously (dashed line in 
Fig. between these power laws but the heuristic spectrum may not be 
correct at intermediate scales since, in particular, Hall MHD turbulence is 
not necessarily anisotropic. However, we will show, in Section [HI that even 
at intermediate scales a moderate anisotropy is expected. The presence of 
a knee in the solar wind spectrum is well attested by in situ measurements 
of magnetic fluctuations (Coroniti et al., 1982; Denskat et al., 1983, Leamon 
et al., 1998; Bale et al., 2005). We will discuss about this point in Section 
ITU1 we will also comment the fact the comparison with observations is not 
direct, in particular, for the high frequency part of the spectrum. One reason 
is that the Taylor hypothesis, usually used at low frequency, is not applicable 
anymore. 
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log(k ± ) 

Figure 2: Sketch of the power law spectrum, at a given k\\, expected for 
incompressible Hall MHD in the regime of wave turbulence. 

7 Nonlocal interactions 
7.1 Introduction 

Nonlocal interactions may play a significant role in turbulence (see e.g. Laval 
et al., 2001). By nonlocal interactions, we mean interactions between well- 
separated scales. For triadic processes, it corresponds to the interaction of 
two highly elongated wavevectors with a third small wavevector. In Hall 
MHD the study of nonlocal interactions and nonlocal transfers are particu- 
larly interesting since the turbulent flow behaves very differently at different 
scales (see Mininni et al., 2005). At the larger scales it follows the MHD 
dynamics whereas at the smallest scales the flow may follow, for example, 
the EMHD dynamics. For simplicity, we will focus our analysis on the mas- 
ter equations (|75p. Therefore we will neglect the helicity effects and keep 
only the kinetic and magnetic energy terms. According to the situation, we 
will consider two different kinds of nonlocal interactions, namely k <C p,q 
or p <C k,q (see Fig. |HJ). In practice, small scale and large scale effects will 
be effective when, respectively, the limits kdi — > +oo and kdi — > will be 
taken. Therefore, we see that large scale effects (small wavenumbers) will be 
driven by Alfven waves (A) whereas small scale effects (large wavenumbers) 
will be driven by whistler (W) or ion cyclotron (C) waves. Note that we 
will not study nonlocal interactions between the same type of waves (AAA, 
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WWW or CCC interactions) which are somewhat more restrictive. Under 
such a condition, three different families of interactions may happen, namely 
WWA, CCA and WCA. We describe below such interactions. 

k « p, q p « k, q 



k 




Figure 3: Nonlocal triadic interactions. 



7.2 WWA interactions 



First we shall consider the strongly nonlocal interactions of two whistler 
waves on one Alfven wave. In other words, it means that the Alfven wave 
will be supported by the wavevector k, with k <C p, q. For such a limit the 
master equations (|7H|) reduce to 



d-, 



\E v {k)\ 



(98) 
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Note that for whistler waves the polarizations are fixed, i.e. A p = s p and 
A q = s q . Further simplifications are possible and one obtains finally 



(E v (k)\ 
\E B (k)S 



(99) 
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E B (q) [E v (k) - E B (k)\ S(Q k , pq ) 5 k , pq dpdq. 

A very interesting result is obtained: at leading order, the Alfven wave is 
only affected by small scale whistler waves when the local (i.e. at a given 
wavevector k) equipartition between the alfvenic kinetic and magnetic ener- 
gies is broken, i.e. when E v (k) ^ E B (k). Additionally, we see that these 
nonlocal effects are the same for the alfvenic kinetic and the magnetic en- 
ergies. This is not so trivial since whistler waves are described only by the 
magnetic energy. We also note that the collisional integral is linearly depen- 
dent on the difference between the large scale kinetic and magnetic energies 
which can be extracted from the integral. These comments may be summa- 
rized by the relation 



h(k) E v (k)-E"(k) 




(100) 



where I\ (k) is a real function that measures the effects of whistler waves. 

The second case that we shall consider is the one where a whistler wave 
is affected by the interaction between an Alfven wave and another whistler 
wave. In other words, it means that the Alfven wave will be supported by 
the wavevector p with p k,q. The same kind of manipulations as before 
lead finally to the following wave kinetic equations 



(101) 




sk — s q q 
, P\\ 
-E v < 



J (sk + s q q) sk\\kp\\ 
5(ttk, pq ) 5k, Pq dp dq. 




K p)-u- [p) 

As expected, we note that at leading order the kinetic energy of whistler 
waves is negligible and only the magnetic energy is relevant. We see that 
the Alfven wave affects the whistler wave only when the local equipartition 
is broken. Thus for any kind of WWA nonlocal interactions, one needs a 
discrepancy from the alfvenic equipartition to have a non trivial dynamics. 



7.3 CCA interactions 

We study here the strongly nonlocal interactions of two ion cyclotron waves 
on one Alfven wave: in other words, it means that the Alfven wave will be 
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supported by the wavevector k, with k <C p,q. For such a limit the wave 
kinetic equations (|75]) write 



a f^ ( k)i 



(102) 
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For ion cyclotron waves the polarizations are A p 
Further simplifications lead to 

aP' (k) l- 



-Sp and A q 
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As before, we see that at leading order the Alfven wave is only affected by 
small scale ion cyclotron waves when the local equipartition between the 
alfvenic kinetic and magnetic energies is broken. Additionally, we see that 
these nonlocal effects are the same for the alfvenic kinetic and magnetic 
energies. This is not so trivial since ion cyclotron waves are described only by 
the kinetic energy. As before, we note that the collisional integral is linearly 
dependent on the difference between the large scale kinetic and magnetic 
energies which can be extracted from the integral. These comments may be 
summarized by the relation 



s p"(k)i 



/ 2 (k) E"(k)-E v (k) 




(104) 



where ^(k) is a real function that measures the effects of ion cyclotron waves. 

The second case that we shall consider is the one where an ion cyclotron 
wave is affected by the interaction between an Alfven wave and another ion 
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cyclotron wave. In other words, it means that the Alfven wave will be sup- 
ported by the wavevector p with p <C k,q. The same kind of manipulations 
as before lead finally to 

E B (k)\ 



(105) 
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As expected, we note that at leading order the magnetic energy of ion cy- 
clotron waves is negligible and only the kinetic energy is relevant. We see that 
the Alfven wave affects the ion cyclotron wave only when its local equiparti- 
tion is broken. Thus for any kind of CCA nonlocal interactions, one needs a 
discrepancy from the alfvenic equipartition to obtain a non trivial dynamics. 



7.4 WCA interactions 



First, we investigate the strongly nonlocal interactions of a whistler wave and 
an ion cyclotron wave on an Alfven wave: in this case the Alfven wave will 
be supported by the wavevector k; the whistler and ion cyclotron waves will 
be associated respectively to p and q, with k p,q. For such a limit the 
wave kinetic equations (|75j) write 



(106) 
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Further simplifications give finally 
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E (q) [E v (k) - E B {k)\ 5(Q k , pq ) 5 k , pq dp dq . 

As before, at leading order the Alfven wave is only affected by small scales 
effects when the local alfvenic equipartition is broken. These small scale 
effects are only due to ion cyclotron waves. The reason of this apparent 
asymmetry is actually due to the initial association between wavevectors 
and waves (see below). We see that these nonlocal effects are the same for 
the alfvenic kinetic and the magnetic energies. Once again, we note that the 
collisional integral is linearly dependent on the difference between the large 
scale kinetic and magnetic energies which can be extracted from the integral. 
Then one can write 



(108) 



where ^(k) is a real function that measures the small scale effects due to the 
ion cyclotron wave. 

Because of the apparent asymmetry mentioned above, we will also con- 
sider the situation where the whistler wave is associated to the wavevector q 
and the ion cyclotron wave to the wavevector p. For such a limit, we obtain 
in the same way as before the following wave kinetic equations 
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The same comments as before hold: at leading order the Alfven wave is 
only affected by small scales effects when the local alfvenic equipartition is 
broken. Here these small scale effects are only due to whistler waves. We 
also see that these nonlocal effects are the same for the alfvenic kinetic and 
the magnetic energies. Once again, we note that the collisional integral 
is linearly dependent on the difference between the large scale kinetic and 
magnetic energies. Then we can write 
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where h(k) is a real function that measures the small scale effects due to the 
whistler wave. 

The next situation that we are going to investigate is the one where an 
ion cyclotron wave is affected by the interaction between an Alfven wave 
and a whistler wave. In other words, it means that the Alfven wave will be 
supported by the wavevector p, the whistler and ion cyclotron waves will 
be associated respectively to q and k, with p <C k,q. The same type of 
manipulations as before lead to 
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As expected, at leading order the magnetic energy of ion cyclotron waves is 
negligible and only the kinetic energy is relevant. Once again, the Alfven 
wave affects the ion cyclotron wave only when the local equipartition is bro- 
ken. 

The last case is the one where a whistler wave is affected by the interaction 
between an Alfven wave and an ion cyclotron wave. In other words, it means 
that the Alfven wave will be supported by the wavevector p, the whistler and 
ion cyclotron waves will be associated respectively to k and q, with p <C k,q. 
Under these conditions and after some simplifications the master equations 
write 
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As expected, at leading order the kinetic energy of whistler waves is negligible 
and only the magnetic energy is relevant. Once again, the Alfven wave affects 
the whistler wave only when the local equipartition is broken. In conclusion, 
for any type of WCA nonlocal interactions, one needs a discrepancy from the 
alfvenic equipartition to obtain a non trivial dynamics. 
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7.5 Discussion 

We have seen that for any kind of nonlocal interactions (WWA, CCA or 
WCA) a non trivial dynamics is present as long as the large scale kinetic 
and magnetic energies are different. This result reveals from the above de- 
tail analysis is new and may explain some solar wind features. Indeed at 
first sight, the solar wind may be seen as a MHD flow where the kinetic and 
magnetic energies are balanced. However, fluctuations around this equipar- 
tition state exist. The present work may explain how these fluctuations are 
maintained by small scales processes. Small scales may also be affected by 
what happens at the largest scales of the system. In this case Alfven waves 
may be seen as a potential source of energy for small scales. We report in 
Appendix O the result of the nonlocal analysis made on WCC interactions. 

8 Source of anisotropy 
8.1 Local interactions 

In previous sections, we have seen three different limits for which anisotropic 
turbulence develops more or less strongly, i.e. the small scale limit, with right 
or left polarization, and the large scale limit of Hall MHD. For intermediate 
wavenumbers the situation is more difficult to analyze. Indeed, the wave 
kinetic equations are in general too complex to give a simple answer. In this 
section, we shall see what happens when strongly local interactions are only 
taken into account. It is a situation where wavevectors have approximately 
the same amplitude, i.e. k « p w q (we consider triangles, k = p + q, 
approximately equilaterals) . We are particularly interested by what happens 
at intermediate scales. 

Our analyze may start from the master equations (|75jl. The hypothesis 
of strongly local interactions does not lead to drastic simplifications; it seems 
therefore more interesting to look at the resonance condition (p>fj) that ap- 
pears in equations (|75j). The first case that we may investigate is the one for 
which waves have the same polarization (left or right). After simplifications, 
one finds the approximate relation 




(113) 



41 



Necessarily two directional wave polarities will be equal which implies that 
one of the parallel wavenumbers will be (approximately) equal to zero. In 
other words, local interactions between the same kind of waves lead to 
anisotropic turbulence where small scales are preferentially generated per- 
pendicular to the external magnetic field. 

The second case that we shall analyze is the one where we have local 
interactions between waves with different polarizations (left and right). Then 
the resonance condition simplifies to 

s p (s p A p +yjl+ A/d 2 p 2 ) - s q (s q A q + y/l + A/ dig 2 ) 

k \\ 

s q (s g A g + sfl + A/d 2 q 2 ) - s(sA + sfl + A/d 2 k 2 ) 

P\\ 

s(sA + + A/d 2 k 2 ) - s p (s p A p + yjl + A/ d 2 p 2 ) 

Several combinations are possible. We note that for the particular case where 
two waves have both the same polarization (left or right) and the same di- 
rectional polarity, one parallel wavenumber will be (approximately) equal to 
zero. Thus, such a situation leads to the same type of anisotropic turbulence 
as before. A similar analysis applied to all the other cases does not lead to 
such a conclusion: an isotropic turbulence may be generated by the nonlinear 
dynamics. 

The wave kinetic equations (J75J) are composed of a sum over the different 
polarities (A- and s-type). In other words, all the cases that we have men- 
tioned above have to be taken into account to understand what happens for 
strongly local interactions. Both the mechanisms leading to anisotropy and 
isotropy are in competition and therefore it is likely that globally a moderate 
anisotropy (less strong than what happens in the limits of large scale or small 
scale Hall MHD) is generated at intermediate wavenumbers. 

8.2 Discussion 

The local interaction analysis gives some hints to understand globally the 
nonlinear dynamics of wave turbulence in Hall MHD. It is also instructive 
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to analyze the opposite limit, namely the limit of strongly nonlocal interac- 
tions. In Section we have introduced the concept of nonlocal interactions 
and discussed about some properties. Information about anisotropy can be 
extracted from the different cases analyzed. More precisely, one needs to 
study the solution of the condition Qk, Pq = that has to be satisfied at the 
resonance. The general tendency that we observe (we will not give a de- 
tail description of every cases) is that most of nonlocal interactions leads to 
anisotropic turbulence where small scales are preferentially generated per- 
pendicular to the external magnetic field. 

Both local and nonlocal interactions contribute to develop anisotropy 
more or less strongly. Additionally, we know that a direct energy cascade is 
predicted in the large scale limit (pure Alfven wave turbulence) as well as 
in the small scale limit (whistler and ion cyclotron wave turbulence) of Hall 
MHD. Thus it is expected that a turbulent anisotropic spectrum generated, 
or developed, at large scale will spread out over wavenumbers, from kdi < 1 
to kdi > 1; extending anisotropy to smaller scales. This scenario may be 
applicable to the solar wind for which various indirect lines of evidence, at 
low and high frequencies, show that waves propagate at large angles to the 
background magnetic field and that the power in fluctuations parallel to the 
background magnetic field is much less than the perpendicular one (Coroniti 
et al., 1982; Leamon et al., 1998). 

9 Wave versus strong turbulence 

9.1 Domain of validity of wave turbulence 

Wave turbulence deals with asymptotic developments which are based on 
a time scale separation, with a transfer time assumed to be much larger 
than the wave period. A consequence of this assumption is that the the- 
ory is not valid uniformly in all of k-space. In this section, we are going 
to evaluate the condition of validity of the present theory. Generally, one 
needs to start with the wave kinetic equations for which exact power law 
solutions exist. In our case, different solutions have been found for different 
limits but in general - for intermediate wavenumbers - no solution has been 
found. For that reason we are going to use a dimensional analysis and, in 
particular, the heuristic spectrum derived in Section |3 This spectrum is 
based on anisotropic arguments, therefore the condition that we are going to 
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write is only valid for anisotropic turbulence and might not be applicable at 
intermediate wavenumbers. 

We remind that the nonlinear time and the wave period are respectively 

TNL ~ Titzi (115) 

and 

11 . . 

Tw ~ 77 = b kw f s ■ (116) 

Note that we have included explicitly the small parameter e in the nonlinear 
time. The transfer time r tr is built on these time scales and writes 

rt r ~rjvi— ■ (117) 
Then the asymptotic condition r tr t w is 

e fc l^A -A)^||£a 
where £ A has a dimension of a velocity (i.e. it is not taken in Fourier space 
as it was introduced in Section [2.5jl . To proceed further we use relations 
an d we take into account the heuristic spectrum (jHSJ); one obtains 
(with B = 1) 

{Zlf ~ (g-gg! fc Z i^(l + fcX)-V4 . (119) 
1+^A 

Then we insert relation (|119|) into equation (|118|) . By using relation (|126|) . 
we finally obtain the inequality 

fey » e 4 / 3 A;i /3 (l + k\d 2 ^l\\ + ef)~ 2/3 • (120) 

Equation ()120|) is the condition of validity of wave turbulence for Hall MHD 
at the level of three- wave interactions. This condition is not in contradiction 
with the anisotropic assumption (k± ^> fen) since the small parameter e is 
present. Thus we can always find a domain of the k-space where the theory 
is valid. We note, in particular, that for whistler and ion cyclotron wave 
turbulence the condition of validity reduces, respectively, to k\\ ^> 

and fey ^> We remind that it is an evaluation based on anisotropic 

assumption and therefore it might not be valid at intermediate wavenumbers. 
In the prohibited region (small fen) other higher-order processes such as four- 
wave interaction have to be taken into account. Basically it is the domain 
where strong turbulence occurs. 
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9.2 The two-dimensional state 

The wave turbulence theory developed in this paper describes the Hall MHD 
dynamics of three-dimensional incompressible waves. Indeed, the condition 
of validity found previously shows that the two-dimensional state, i.e. modes 
with fey = 0, is not described by the wave kinetic equations. Actually, if we 
look at equations (|56|) . we see that the nonlinear transfers decrease (linearly) 
with k\\. For the forbidden value k\\ = the transfer is exactly null. We see 
therefore that the two-dimensional - slow - modes decouple from the three- 
dimensional Hall MHD waves. Such decoupling is not new; it is found in a 
variety of problems like for internal gravity waves (see e.g. Phillips, 1969), 
rotating stratified flows (see e.g. Bartello, 1995) or pure rotating flows (see 
e.g. Smith and Waleffe, 1999). 

As we have already noticed in Section 14.51 the MHD limit is somewhat 
singular for Hall MHD. A signature of this singular behavior is the appear- 
ance of principal value terms in the wave kinetic equations. The consequence 
is that there is an important difference between Hall MHD and MHD: in the 
latter case the two-dimensional modes drive the three-dimensional dynamics. 
Therefore, in this particular case, we need information from the forbidden 
region to describe the dynamics in the domain of applicability of wave tur- 
bulence. Note that even in this case the wave turbulence theory at the level 
of three- wave interactions has been shown to be applicable (Galtier et al., 
2000; Nazarenko et al., 2001, Bhattacharjee and Ng, 2001). For Hall MHD 
in general, higher order processes like four-wave interactions may lead to a 
coupling between the two-dimensional state and three-dimensional modes. 



9.3 Strong Hall MHD turbulence 

We shall derive now the equivalent heuristic prediction as in Section |H1 but 
for strong turbulence. We would like to compare, at least at the level of the 
phenomenology, the wave and strong turbulence regimes. 

We will assume that a moderate external magnetic field B is applied 
and that it still affects the Hall MHD dynamics. In other words we will 
distinguish the wavenumber k± from k\\ , with k± S> k\\ , and we suppose that 
the nonlinear time is of the order of the wave period. Then one has 

n~ *~ g ^*^y (121) 

T~tr T NL 
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The equivalent manipulations as in Section lead to the prediction 

E{k ± , k\\) ~ n 2 / 3 fcj^jfejj-^l + 5 As k ± di)- 2 / 3 , (122) 

where 5\ s is introduced to distinguish the electron MHD case (5 SS = 1) from 
the ion MHD case (5- ss = 0). The heuristic prediction gives the energy spec- 
trum scaling laws for the three different limits of standard, electron and ion 
MHD as well as for intermediate wavenumbers. Contrary to the wave tur- 
bulence regime the prediction shows that EMHD and IMHD follow different 
power laws. If the external magnetic field is small enough one may ignore 
anisotropy and write k± ~ k\\ ~ k\ it gives 

E(k) ~ n 2 / 3 k~ b /\l + 6M)- 2 / 3 , (123) 

for the one- dimensional isotropic spectrum. It is interesting to note that the 
EMHD prediction is the one found and observed in direct numerical simu- 
lations by Biskamp et al. (1999). Note that the Kolmogorov scaling found 
at large scales is an unavoidable result since we have used in fact the same 
phenomenology as Kolmogorov. Note finally that the scaling law predic- 
tions made for strong and wave turbulence in the presence of anisotropy may 
be described in the context of the generalized critical balance proposed by 
Galtier et al. (2005). In that context the solutions found here would be two 
particular cases of a familly of solutions. The result obtained here has to 
be contrasted with those obtained by Krishan and Mahajan (2004,2005) for 
strong isotropic turbulence. In these papers, the authors give several predic- 
tions separately for the magnetic and kinetic energy spectra. In particular, 
they conclude that the whistler-type prediction is very different from the ob- 
servations and evoke a possible stronger damping for whistler waves. In the 
present paper, it is thought that the damping (e.g. due to the resonance) 
affects mainly the ion cyclotron waves (see also Goldstein et al., 1994). 

10 Conclusion 
10.1 Summary 

In this paper we have investigated the steepening of the magnetic fluctua- 
tion power law spectra observed in the solar wind for frequencies higher than 
0.5 Hz. We have shown that the high frequency part of the spectrum may 
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be attributed to dispersive nonlinear processes present in the incompressible 
Hall MHD equations rather than pure dissipation. In that context, we have 
introduced the generalized Elsasser variables adapted to Hall MHD and de- 
veloped a wave turbulence formalism based on a complex helical decomposi- 
tion. Then, we have derived the wave kinetic equations for three dimensional 
incompressible Hall MHD turbulence at the level of three-wave interactions, 
in the presence of a strong external magnetic field, and in the most general 
case, i.e. with global and magnetic helicities. In such a regime, Hall MHD 
turbulence is mainly anisotropic with small scales preferentially generated 
perpendicular to the external magnetic field. In particular, strong anisotropy 
is developed at small scales (ion cyclotron/IMHD and whistler /EMHD wave 
turbulence) and large scales (pure Alfven/MHD wave turbulence). For inter- 
mediate wavenumbers, i.e. kdi ~ 1 only a moderate anisotropy may happen. 
We have found two different power law energy spectra that are exact solutions 
of the wave kinetic equations, with a steeper power law at smaller scales. To 
illustrate this finding, we have developed an anisotropic phenomenology that 
describes continuously the different scaling laws for the energy spectrum: the 
existence of a double inertial range is recovered with a knee at intermediate 
scales where the Hall term becomes (sub-) dominant. We have also shown 
that the large scale limit of Hall MHD, i.e. the standard MHD, is singular 
and leads to the appearance of principal value terms in the wave kinetic equa- 
tions. We have analyzed the nonlocal interactions between Alfven, whistler 
and ion cyclotron waves and we have shown that, at leading order, a non 
trivial dynamics is present only when a discrepancy from the equipartition 
between the large scale kinetic and magnetic energies of Alfven waves hap- 
pens. Finally, we have given the conditions of validity of wave turbulence in 
Hall MHD and compared the wave and strong turbulence scaling laws. 

10.2 Discussion 

Solar wind magnetic fluctuation power law spectra are commonly observed 
to have an index around —5/3 at frequencies lower than 0.5Hz and a steeper 
power law behavior at higher frequencies with a spectral index on average 
around —3 (Coroniti et al., 1982; Denskat et al., 1983, Leamon et al., 1998). 
The latest analysis made from Cluster spacecraft data has - however - re- 
ported a spectral index about —2.12 (Bale et al., 2005). These small scale 
fluctuations are also characterized by a bias of the polarization suggesting 
that they are likely to be right-hand polarized, outward propagating waves 
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(Goldstein et al., 1994). For these reasons, it is claimed that Alfven - left 
circularly polarized - fluctuations may be suppressed by proton cyclotron 
damping and that the high frequency power law spectra are likely to consist 
of whistler waves (see e.g. Stawicki et al., 2001). Additionally, it is often 
observed that the high frequency fluctuations of the magnetic field are much 
smaller than the background magnetic field (see e.g. Stawicki et al., 2001). 
Therefore, a wave turbulence approach seems to be well adapted to this 
problem. It is well known that whistler waves are solution to the Hall MHD 
equations. We have seen in the regime of wave turbulence that when the 
Hall term is added to the MHD equations a nontrivial dynamics happens 
at small scales (i.e. when kdi > 1) with, in particular, a steepening of the 
power law energy spectra. This behavior is similar to what is observed in 
the solar wind and may be seen as an indication that the steepening of the 
solar wind magnetic fluctuation power law spectra is mainly due to nonlinear 
processes rather than pure dissipation. Under this new interpretation, the 
resistive dissipation range of frequencies may be moved to frequencies higher 
than the electron cyclotron frequency. 

The Taylor hypothesis, that allows to connect directly a frequency to 
a wavenumber, is widely used to interpret the single spacecraft solar wind 
data. It is thought that this approximation is well adapted to analyze, in 
particular, the low frequency part of the magnetic and velocity fluctuations. 
The subsequent interpretation is mainly relevant for isotropic media. How- 
ever, there are evidences that anisotropy is present in the solar wind at high 
and low frequencies (see e.g. Belcher and Davis, 1971). From a theoretical 
point of view, it is well known that the presence of a strong magnetic field 
influences the MHD turbulent flows (see e.g. Pouquet, 1978; Montgomery 
and Turner, 1981; Shebalin et al. 1983; Goldreich and Sridhar, 1995; Ng 
and Bhattacharjee, 1996-1997; Galtier et al., 2000-2005). The main effect is 
that MHD turbulence becomes mainly bidimensional with a nonlinear trans- 
fer essentially perpendicular to its direction. Direct numerical simulations 
of 2|D compressible Hall MHD for high and low beta (3 plasma (Ghosh and 
Goldstein, 1997) have also displayed such an anisotropic property when a 
strong magnetic field is present. As we have seen, the model that we propose 
here is also able to exhibit such an anisotropy. The predictions that we have 
made are for anisotropic turbulence, i.e. a situation where we distinguish the 
wavenumber k± from k\\ . Therefore, a proper comparison with observational 
data will be possible only when a three dimensional energy spectrum will be 
accessible. It is important to note that average effects may alter significantly 
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the power law scaling. An illustration of such effects may be given from 
the recent heuristic MHD predictions made by [39j who have generalized 
the concept of critical balance. For a medium where a strong or moderate 
magnetic field Boen is present, they predict the anisotropic energy spectrum 
E(k±, k\\) ~ k2 a k^, with 3a + 2(3 = 7. This model is able to describe both 
the strong and wave turbulence regimes as well as the transition between 
them; it also satisfies the critical balance relationship ku oc k^ 3 . It is very 
interesting to note that if one averages the heuristic spectrum to obtain the 
one dimensional counterpart E(k±), one obtains systematically a scaling law 
in kj_ for any family of solutions (a, (3). Actually, this remark might ex- 
plain the apparent contradiction between, on the one hand, the presence of 
anisotropy in the solar wind and, on the other hand, the Kolmogorov (—5/3) 
scaling law for the low frequency energy spectrum. Of course, other explana- 
tions of the 5/3-spectrum are possible; an example is given in Oughton and 
Matthaeus (2005). The direct comparison between theoretical predictions 
and in situ measurements of three dimensional energy spectra is particu- 
larly crucial for the high frequency part of the power spectrum for which 
the usual Taylor hypothesis that allows to connect directly the frequency to 
the wavenumber is not applicable anymore. Efforts are currently made with 
Cluster spacecraft data from which it is possible to extract the three dimen- 
sional magnetic turbulent spectra of the magnetosheath thanks to multipoint 
measurements and a k-filtering technique (see e.g. Sahraoui et al., 2004). As 
explained in j^j, Cluster spacecraft may exit from the terrestrial magneto- 
sphere to make solar wind measurements. The application of the k-filtering 
technique to the high frequency part of the solar wind magnetic fluctuations 
seems then possible. It may lead, for the first time, to a direct and rigorous 
comparison with a model prediction of solar wind. Note that using multi- 
spacecraft analysis, Matthaeus et al. (2005) have recently investigated some 
spatial correlations from two-point measurements. 

The model proposed in this paper is able to recover some solar wind prop- 
erties but several aspects have not been discussed. For example, it would be 
interesting to investigate the effects of asymmetry, i. e. the fact that outward 
propagating waves dominate in the solar wind. In the framework of incom- 
pressible MHD, we know that asymmetry changes the index of the power 
law spectra. Similarly, the indices found here may be affected by asym- 
metry; this question will be tackled in the future. In this paper, we have 
found that a turbulent state made of ion cyclotron waves may exist around 
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a fraction of the ion cyclotron frequency Ud, namely for a frequency around 
oj r = (k\\/k) ijJ C i- The presence of a resonance at a frequency uir lower than 
ojd is rarely discussed in the litterature (see White et al., 2002) since the 
analysis is generally focused on parallel propagations (Tcii = k) for which the 
resonance frequency is exactly the ion cyclotron frequency As we have seen, 
in the presence of a strong external magnetic field, Hall MHD turbulence be- 
comes mainly bidimensional with an energy spectrum mainly spreaded out 
over the wavenumbers k±. In this case, the resonance frequency may appear 
at a small fraction of u C i (this fraction being even smaller for heavier mass 
ions). These remarks may be crucial to understand the solar coronal heating 
problem in which the coronal temperature is far beyond what one can pre- 
dict by the resistive MHD approximation: although the Hall MHD model is 
a fluid model that does not describe the resonance between waves and partic- 
ules and therefore the particule heating, it offers the possibility to evaluate 
the rate of particule heating by assuming that the turbulent energy of ion 
cyclotron waves are mainly transfer into heating. This point is currently in- 
vestigated and will be reported elsewhere. In the context of the solar wind, 
the measure of the position of the knee in the magnetic fluctuation power 
law spectrum may be seen as a proxy to measure the solar wind anisotropy. 

The Hall effect is relevant in many astrophysical problems to understand, 
for example, the presence of instabilities in protostellar disks (Balbus and 
Terquem, 2001), the magnetic field evolution in neutron star crusts (Gol- 
dreich and Reisenegger, 1992; Cumming et al, 2004) or impulsive magnetic 
reconnection (see e.g. Bhattacharjee, 2004). Small scale turbulence is also 
a key issue in a number of problems from the interstellar medium (see the 
review made by an d, jHH]), to solar physics (see e.g. , |85j), magneto- 
spheric physics (see e.g. the recent paper by jlH]) and laboratory devices such 
as tokamaks (see e.g. Wild et al., 1981; Taylor, 1993). Therefore, it is likely 
that the present model will interest several other (astrophysical) problems. 
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A Some useful relationships 

It is convenient to note the following identities: 

l-a 2 = kAdi&, (124) 

£ur = -i, (125) 



a - = -«v 4 + ^ ' ( 126 ) 
a+a s = -^A, (127) 



A l^—^-l I -Z^-j 1=4*11. (128) 



A 



B Derivation of the wave kinetic equations 

The starting point of the derivation of the wave kinetic equations for incom- 
pressible Hall MHD is the fundamental equation (|4*7|) . We write successively 
equations for the second and third-order moments, 

dt(a' A a&) = (129) 

£a 9 - AApAq 
j-f / E a 2 7_ ' ,- P s M (a£a£a£> e ^.*% >Jk cfedq 

+ 

J a p ,a 9 ?a' — ?A' ^ i^y 

Sp.Sq 

and 

<9 t (a A a A ',a A ",) = (130) 

t«« _ t«f AA P A, 

^U>i J A„.A„ tA CA ft X) 



^tdi J A „ A _ C A / C A / K T) Q 



+ 

•j> A'ApA, 



"l Ad, A 
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+ 



e 



"2 SA 



t s 9 dr s p 

SA„ SA„ 



SA" SA" 



M 



-k" pq 



(«A «A «A«A') e iUpq ' kl,t <W" dp dq 



We shall write asymptotic closure (Newell et al., 2001) for our system. For 
that, we basically need to write the fourth-order moment in terms of a sum of 
the fourth-order cumulant plus products of second order ones. The asymp- 
totic closure depends on two ingredients: the first is the degree to which 
the linear waves interact to randomize phases; the second relies on the fact 
that the nonlinear regeneration of the third-order moment by the fourth- 
order moment in equation (J 130)) depends more on the product of the second 
order moments than it does on the fourth order cumulant. The fourth- 
order moment decomposes into the sum of three products of second-order 
moments, and a fourth-order cumulant. The latter does not contribute to 
secular behavior, and among the other products one is absent because of 
the homogeneity assumption. If we use the symmetric relations (I49j) - ()52j) 
and perform wavevector integrations, summations over polarities and time 
integration, then equation (jl30|) becomes: 



Ad, 



^{Qkk'k") &kk'k" 



(131) 



s 2 
A 



SA" ~ SA' 



AA'A" \ * 

t a" \ 



M ss ' s ' 

h h' h" 



A'AA" - 

' 6 »" 



where 




ts' _ ts" / 
_j_ SA' SA" M Ss - s 



AA"A' \ * 



k k" k> 



+ 



A'A"A 



e in kk'k"t'dt' 




iVthu' 



_ £a" I jy- a > S " s 

a''-G s 'v k>k«k) 
+ 

ts' _ts / A"AA' - *' 

+ „ „ \M s " ss ' 



Qa' Qa" 



Qa 9A" 



Qa' ll } 



(132) 



kk'k" 



53 



The introduction of symmetric relations (|49|) -(|52 |) into (|131|) allows us to 
simplify further the previous equation; one obtains: 



(a S A a S A ,a S A ,,) — A(Q kk > k 't) Skk'k'' 



AA'A" 



JY£ s s' s" 

h y y 

IX IX fx 



(133) 



C s2 ZA" ~ S"A' „s' s" j_ ts' 2 S~A ~ SA" „s s" , C s" 2 ^A' ~ >A „s *> 
SA Ts _ ==s "A' <?A" "+" SA' TTv^A VA" "+* SA" TZ^^A HA> 

SA SA SA' SA' SA" SA" 

We insert expression (|133|) into equation (|129|) : it leads to: 
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—kpq 
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£Sq £Sp 

2 SA„ _ SA D 



SA D £ sp f -s p «A' VAg + SA, c s q c -sg HA' HA P + SA' ts' _ t-s> HA V HA 
SAp — SAp SA„ — SA„ SA' SA' 



dp cfq . 



The long-time behavior of the wave kinetic equation ()134|) is given by the 
Riemman-Lebesgue Lemma which tells us that, for t +00, we have 



-ixt 



A(x) = A(-x) -> tt8(x) - iV(l/x) 



(135) 



where V is the principal value of the integral. The two terms of equation 
(I134j) are complex conjugated therefore if in the second term we replace 
the dummy integration variables p, q, by — p, — q, we can simplify further 
equation ()134|) since, in particular, principal value terms compensate exactly. 
Finally, we obtain the wave kinetic equations for incompressible Hall MHD: 



(136) 
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where 



AApAq 2 



j\/f SS P S 9 

—kpq 



The last step that we have to follow to obtain the same expression as (|56j) is 
to include the resonance relations (jHH) into the previous equations. 

C Pseudo-dispersive MHD waves and com- 
plex helical decomposition 

This section is devoted to the derivation of the wave kinetic equations for 
pure alfvenic turbulence (kdi = 0). These equations were already derived by 
Galtier et al. (2000) but here we use the complex helicity decomposition. 
We will see that some differences appear in the kinematics that renders the 
MHD description somewhat singular (principal value terms appear) by op- 
position to the dispersive Hall MHD description. We start from the standard 
incompressible and inviscid MHD equations, 



V ■ V = 0. 



— + V • V V = -VP* + B V B 



«9V 

~dt 



dB 

~dt 



+ V VB = B - VV 
V-B = 0. 



(137) 

(138) 

(139) 
(140) 



The same notation as before is used. We introduce the fluctuating fields 
B(x) = B e\\ + eb(x), V(x) = ev(x) and we Fourier transform the MHD 
equations. One obtains: 



d t z\ - isk^zX = -e {z~ s • V z s - VP*}! 



141^ 
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k-z s k = 0, 



(142) 



where we use the standard Elsasser variables z s = v + sb. The linear solution 
(e = 0) corresponds to linearly polarized Alfven waves for which uj^ = B kn. 
We introduce the complex helicity basis and define 

< = E 2i(k) e ^ K = Y. z l^ kt K ■ (143) 

A A 

We see that this decomposition is not natural for linear polarized Alfven 
waves since for a given direction of propagation (a given s) we have two 
contributions for each value of A. We nevertheless use this decomposition to 
show the compatibility with the Hall MHD turbulence description. After the 
substitution of (11331) into (IT4T1) - (IT4"21) . we obtain 

d t Zi = ~ f £ Z A ;^(k-h^)(h^.h k A )e-^ + ^)' ( 5 w , fc rfp C /q. 

Z J A p ,A q 

(144) 

Note that we do not have a summation over the directional polarity s p and s q . 
The physical reason is that the nonlinear coupling in incompressible MHD 
involves only Alfven waves propagating in opposite directions. Thus the 
information about the direction of propagation is already taken into account 
in equation (|141|) . We use the local decomposition and find after some algebra 

d t Z A = jf E ( kA + <? A <? - pA P ) Z H Z k e ~ 2lSB(mt W dp dq , 

4 J A P ,A g 

(145) 

where 

N^ q A " = e *(A* fc +A P P +A g * ? ) A Ap Ag (AA . + App + Aqq) (146) 

We note that the matrix iV possesses the following properties (* denotes the 
complex conjugate), 

* = X-tf-t* = At , (147) 
N^f p , (148) 

KlT » (1 49 ) 
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Equation (|145|) is the fundamental equation that describes the slow evolu- 
tion of the Alfven wave amplitudes due to the nonlinear terms of the incom- 
pressible MHD equations. Note that (i) we have already used the resonance 
condition to simplify the coefficient in the complex exponential function; (ii) 
a comparison with the large scale limit of equation (jTFj) is possible if we sum 
over the directional polarities. We follow the same steps as in Appendix [B] 
we write successively equations for the second and third-order moments, 

d t (ZlZi) = (151) 
- f E ( kA + ?A, - pAp) N*t P p A q q (Zl°ZiZi)e- 2isB °nt Spq , k dp dq , 

Ap,A 9 

+ 

7 / £ (*'A' + qA q - PK) N-££ (ZgZ^Z' K )e- MB «n* 5 pq , k , dp dq , 
and 

d t {Z s A Z(,Zf„) = (152) 
- ( £ (kA + qA q - P A P ) N A t;f (Z^ZtZiZ^e- 2 ^^ 5pq , k dp dq , 

4J A P ,A 9 

+ 

7 / £ (k'A' + qA q - P A P ) N*£fi (Z^ZiZiZ%,)e- 2is ' B °nt 5pq>k , dp dq , 

A P ,A„ 

+ 

7 / £ (W' + gA 9 - P A p )<'^ 

A P ,A 9 

We define the density tensor g^,(k) for an homogeneous turbulence, 

(Z s A (k) Z&M) = q s A Uk) 5(k + k') 6 SS , . (153) 

The presence of the delta 6 SS ' means that correlations with opposite polarities 
have no long-time influence in the wave turbulence regime; the second delta 
distribution <5(k + k') is the consequence of the homogeneity assumption. We 
note that the kinematics does not impose any condition about the polariza- 
tion A. The reason is that A does not appear in the Alfven wave frequency. 
This remark shows a fundamental difference with the Hall MHD case. As we 
will see it is the reason why principal value terms appear in MHD but not 
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in Hall MHD. After the same kind of manipulations as in Appendix [0 the 
third-order moment equation (|152j) becomes: 

(Z' A Ztzg,) = (154) 

\ <w{ E ( kA + k "K - k ' A v) { N Tk'k")* A{2sBok ^ q ^ A ' ql > 
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We insert expression (jl54j) into (jl51|) . We note that only the third and the 
fifth terms of expression (jl54j) will contribute. One obtains 

«A'(k)= (155) 

U f E (*A + qA q - pA p ) (qA q + kA q - pA p ) N AApAq (-N* 
1,3 J A P ,A 9 —Kpq \ 



AqApAgY 

—kpq I 



Ap,A„ 



A{-2sB m ) q A y qf AI 8 pq>k dp dq 



+ 
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/ E (*A + ,A, - P A P ) (ktt + fA, - P A p) ( N ^) 

Ap.Aq 

A(-2sB oPll ) q^ qf qAq 5 pq , k dp dq 
+ 



A'A„A„ / ..AgApAg 



/ £ <*A< + ,A, - P A P ) (q A, + *A, - P A p) N ' • f -AT 



A(-2sB„pn) (,''jf A 4„<iprfq 



+ 



/ J2(kA' + qA q -pA p )(kA + qA q -pA p )N 



A'A„A„ f N AA p A q s 
kpq 



A(-2sB oPll ) g A ; Ap S q s K s qAg 5 kpq dpdq}. 
We change signs for wavevectors p and q in the last two terms and obtain: 



^Cv(k) = ^/ £ (156) 



{ (kA+qA q - P A p ) (qA q +kA q - P A p ) (^-N A ^J A(-2sB oP]l ) q^q^ 

+ 

(kA+qA q -pA p ) (kA'+qA q -pA p ) (n^J A(-2sB oP]] ) q^ 

+ 

(kA'+qA q -pA p ) (qA q +kA q -pA p ) [~N^J A(2sB oPll ) qjtf q- A 

+ 

(kA'+qA g -pA p ) (kA+qA q -pA p ) N A _ A k P ^ (N A ^J A(2sB oPll )q A X<lA q A q } 

5 mk dp dq . 
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The long-time behavior is then given by the Riemman-Lebesgue Lemma: 

A(±2sB p l{ ) -> n5(2sB QPll ) T iP(l/2sB oP \\) . (157) 

We see that principal value terms will appear in the long-time limit. The 
reason is that the polarization A and A' in the density tensor q AA , (k) are not 
the same in general. Therefore we lose the symmetry between terms that we 
had in the Hall MHD case (see Appendix [Bj. 



D Simplified MHD wave kinetic equations 

In in section, we continue the analysis made in Appendix IU1 when only terms 
symmetric in A are retained, i.e. terms like q AA . Then expression 
simplifies to 



dtq S A(k) 



16 



E 



(158) 



{(kA + qA q - P A p ) 2 N AA u pAc! 

—kpq 



~ NA -kpqj A (~ 2SB ^ Ca7 ?AA 



(kA + qA q - pA p ) 2 N 



AA n A 



p ls -q 



N 



{kA + qA q -pA P YN AKpK « 

—kpq 



+ 

AA p A q 
—kpq \ —kpq 

+ 

AA„A 



A(-2sB p {l ) Qa, 



—kpq 

+ 

.AA„A 



A(2sB oPll ) qXX q s A 



{kA + qA q -pA p fN A ^ q [N^ q 

—kpq V —kpq 

<W dp dq . 



A(2sB oPll ) q^X QZa,} 



The symmetry between terms is recovered; further simplifications lead to: 



16 



E 



^AApAg 

—kpq 



(kA + qA q - pA p 



(159) 



(A(-2sB oPll ) + A(2sB p l{ )) 9 ££ (q s A s qAq - q s A \) <W dp dq, . 
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The long-time behavior is given by the Riemman-Lebesgue Lemma; one finds: 




dtq s A(k) = (160) 
(kA + pA p + qA q ) 2 (kA + qA q - pA p ) 2 

QaX (<A 9 - Q'A) 5 (P\\) 5 k, Pq dp dq . 

Equations ()160|) are the wave kinetic equations for incompressible MHD tur- 
bulence at the level of three- wave interactions when helicities are absent and 
when equality between shear- and pseudo-Alfven wave energies is assumed 
(see Section 14. 5|) . We find the same equations as (J72j) since we have the 
relation 5^(k) = 4g^(k). 



E Kinetic equations for the energies E and 



E. 



d 



We introduce the expression (|71j) into (}56j) and we consider a state of zero 
helicities (H m (k) = and Hoik) = 0). The kinetic equations for the energies 
E and Ed are 

' EQc) 1 



(161) 



7r e 



32 dfBt 



E 



' sin ip k 
k 



2 (Ak + a p p + A q q) 2 i-er exr e 



,._.2 / .- s „2. 



:i+a s2 )(i+c p2 )( i +c q2 ) 



'e:;-ei;y rer+i' 



£(p) - £(k) + 



g +1 

\ ^An ^ / 



SA 
£<i(p) 



£(q) + 



+ i 



^(q) 



^d(k) 



5(^fe,p 3 ) h, Pq dp dq. 
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F Compatibility with Galtier et al. (2000) 

We start from equation (26) of Galtier et al. (2000). We assume that helici- 
ties are absent (/ s (k) = 0) and that shear- and pseudo-Alfven wave energies 
are identical. We introduce the following notations: 

k 2 ± ^ s (k) = e s (k)/2, (162) 

and 

k 2 ± k 2 $ s (k) = e s (k)/2. (163) 

Then equation (26) writes: 

d t e s (k)= (164) 



n e 



2 r (k x k) 2 kin 2 , (k x k) 2 kin 2 . 



45 7 LV L^/c 2 L\k 2 v 7 v £; 2 L 2 fc 2 L 2 

+ 

kf,n 2 cos 2 ib k kf,K 2 cos 2 iIjl 

(1 - " Llk2 ) e s (L) - (1 - " ^ ) e*(k) } 

{ (k x k) 2 - k 2 K 2 ± ) 6 + k 2 e~ s (n)} <5(«||) 5 ktKL dK,dL . 
k ± 

We remind that the angle ipk refers to the angle opposite to the 3D vector 
k in the triangle defined by k = L + k. (In Galtier et al. (2000), angles are 
introduced in reference to 2D wavevectors.) Some manipulations lead to: 

d t e s (k) = (165) 



7r e 



2 fr< k\ sin 2 il) L k 2 n\ kuK 2 ± cos 2 ^ k 



{(2 - ^TY^ + 7^ " " „ , 2 ) eS ( L ) 



ABj LV L 2 : L 2 k 2 L 2 k 



(k 2 Klsm 2 jj L k 2 K 2 ± k\K\ cos 2 ipL 

rvj^lj rvj_±j fX'j^-L/ 

k 2 sin 2 ip L e~ s (K,)5(K,\\) 5 kjKL d,KdL , 
that can be written as: 

<9 t e s (k) = (166) 

2 

/{(l + cos 2 ^ K ) L 2 sin 2 V fc e" s (/t) (e s (L) - e s (k)) 6(k\\) 5 KKL dndL . 
4-Do J 
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Equations (|166|) are the wave kinetic equations for incompressible MHD tur- 
bulence when helicities are absent and when equality between shear- and 
pseudo-Alfven wave energies is assumed. The addition over the index s will 
give the wave kinetic equations for the total energy. 



G Nonlocal WCC interactions 



We shall consider the strongly nonlocal interactions of two ion cyclotron 
waves at small scales on one whistler wave at large scales. In other words, 
it means that the whistler wave will be supported by the wavevector k, with 
k <C p,q. This type of interaction is the most interesting since smaller scales 
may be reached more easily by ion cyclotron waves. For such a limit, the 
master equations (|75|) reduce to 



(167) 




(s p p + s q q) 



2 pq 2 k 




£ y (q) 



(p) - E (k)J 5(tt km ) 5 km dp dq . 

As expected, at leading order the whistler wave is described by the magnetic 
energy. Contrary to the other cases studied before, a non trivial dynamics 
happens as long as a discrepancy exists between the (kinetic) energy of the 
ion cyclotron waves and the (magnetic) energy of the whistler wave. 
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